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ASYMPTOTIC  THEORY  OF  TURBULENT  WALL  JETS 


R.E.  Melnik* 

A.  Rubel** 
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Grumman  Aerospace  Corporation 
Bethpage,  New  York  11714 


Abstract 


This  paper  presents  a  systematic  analysis  of 
two-dimensional  and  radial  turbulent  wall  jets  using 
the  method  of  matched  asymptotic  expansions.  The 
asymptotic  solution  is  carried  out  in  two  stages. 

The  first  is  based  on  a  two  parameter  expansion  of 
the  full  Reynolds  averaged  equations  with  a  k-c 
model  of  turbulence  quantities.  One  of  the  small 
parameters ,  y ,  is  related  to  the  nondimensional 
friction  velocity,  uT,  defined  by  the  surface  shear 
stress.  It  is  basically  a  Reynolds  number  param¬ 
eter,  y  =  0(ln  Re)~*.  that  primarily  controls  the 
shear  stress  effects  induced  by  the  wall.  The 
other  small  parameter,  a,  is  related  to  the  modeling 
constants  arising  in  the  chosen  turbulence  closure. 
In  the  present  k- 1  model  analysis  «  is  identified 
with  the  constant,  cp  ,  appearing  in  the  eddy  vis¬ 
cosity  formula  for  the  Reynolds  shear  stress  and  is 
a  measure  of  the  turbulence  levels  in  the  outer  free 
jet  part  of  the  flow.  The  »  expansion  reduces  the 
problem  to  a  classical  boundary  layer  formulation  to 
lowest  order.  The  expansion  for  y->-0  leads  to  a 
four  layer  description  of  the  wall  jet.  The  outer 
layer  is  closely  related  to  a  free  jet  flow  while  the 
innermost  layer  is  a  classical  law  of  the  wall  region. 
Two  additional  intermediate  layers  are  needed  to 
effect  the  matching  of  the  outer  and  inner  layers 
and  to  complete  the  solution.  Leading  order  solu¬ 
tions  for  each  layer  are  presented  and  the  compos¬ 
ite  flow  field  result  is  compared  with  an  existing 
numerical  solution  for  the  wall  jet. 

Nomenclature 


a  Constant  Defined  by  Eq.  (45b) 

A(x)  Function  Appearing  in  Inner  Limit  of 
Outer  Solution 

As  Constant  Appearing  in  Inner  Limit  of 
Outer  Similarity  Solution 
b  Constant  Given  by  b  =  1  -  (3/2)a 

b*  Reference  Length  in  Initial  Jet 
B(x)  Function  Appearing  in  Inner  Limit  of 
Outer  Solution 

B  Function  Appearing  in  Inner  Limit  of 

s  Outer  Similarity  Solution 

Cf  Skin  Friction  Coefficient  Based  on  Initial 

Jet  Conditions 

Cp  Skin  Friction  Coefficient  Based  on 

uw(x) 

cw(x)  Function  Defined  by  Eq.  40 

<CW  = Uw  duw/dx) 

c(x)  Function  Appearing  in  Transition  Layer 
Solution,  Defined  by  Eq .  (80a) 
ca  ce2  Turbulence  Model  Constants  (c  .  —  1.44. 
c  =  1.92) 

L  L 

c.  Function  of  P/L  and  f  Defining  Turbulent 

u  Viscosity 
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So  Turbulent  Model  Constant  (cj;o  =  0.09) 

C]  Momentum  Flux  Constant  Defined  by  Eq. 

(32) 

C+  Constant  in  Law  of  the  Wall,  Eq.  (103) 
d  Constant  Defined  by  Eq.  (100) 

e  Constant  Defined  by  Eq.  (101) 

E(;)  Similarity  Function  for  Dissipation, 

Transition  Layer  Eq.  (84) 

E  .  Constant  in  Similarity  Solution  for  Dissipa¬ 
tion,  Transition  Layer 

f  Wall  Damping  Function  Defined  by  Eq.  (3) 

F(n)  Streamfunction  in  Similarity  Solution, 

Outer  Layer,  Eq.  (25)  -  (26) 
g  Normalized  Reynolds  Shear  Stress  Defined 

by  Eq.  (10c) 

G(n)  Similarity  Function  for  Turbulent  Energy  - 
Outer  Layer,  Eq.  (27) 

G  Function  Defining  Turbulent  Viscosity  in 

Transition  Layer,  Eq.  (86) 

G„  Function  Defining  Turbulent  Viscosity 

Eq.  (2) 

GMW  ,Gy„  Limit  Values  of  Turbulent  Viscosity 

Function  Defined  by  Eqs.  (43)  and  (47) 

H(n)  Similarity  Function  for  Dissipation  -  Outer 
Layer,  Eq.  (29) 

j  Index  Equal  to  0  for  Plune  Jet,  Equal  to 

1  for  Radial  Jet 

J  Momentum  Flux  Defined  by  Eq.  (24) 

k  Normalized  Turbulent  Energy,  Eq.  (lOd) 

K(n),K(  ,)  Similarity  Functions  for  Turbulent 

Energy  -  Outer  Layer  Eq.  (28)  and  Transi¬ 
tion  Layer  Eq.  (83) 

P  Turbulent  Production,  vt(3  8/3y)^ 

r  Nondimensional  Coordinate  along  the  Wall 

(r  =  r  */b*) 

Rej  Jet  Reynolds  Number  (uj*b*/v*) 

u  Velocity  in  the  x.  r  Direction  -  Outer 

Solution  Eq.  (17) 

uw  Velocity  at  the  Wall  (y  =  0)-  Outer  Layer 

Solution  _  _ 

uT  Friction  Velocity  y\w  =  J  Cf/2 

U  Velocity  Component  in  x,  r  Direction, 

Eq.  (10a) 

Um  Maximum  Velocity  in  Streamwise  Velocity 

Profile 

v  Velocity  in  y  Direction  -  Outer  Solution, 

Eq.  (18) 

V  Nondimensional  Velocity  Component  in  y 

Direction,  Eq.  (10b) 

x  Stretched  Coordinate  in  r  Direction,  Eq.  (11) 

y  Nondimension  Coordinate  Normal  to 

Plate,  Eq .  (11) 

ym  Coordinate  of  Velocity  Maximum 

t  Small  Parameter  Equal  to  cM0 

,  (K)  Function  Appearing  in  Transition  Layer 
Similarity  Solution  Eq.  (80) 

•w  Value  of  6  at  the  Wall,  y  =  0,  Eq.  (107) 

V  Small  Parameter,  y  -  u1>r/  i(1+l)/4 

'  •  :  m  Y  Coordinate  Value  of  Wall  Jet  Edge. 

:  1/2  Maximum  Velocity,  and  Half  Velocity  Point 
j(x).  Functions  Appearing  Diffusion  Layer 

Solutions.  Eqs.  (63)  (64) 

Normalized  Turbulent  Dissipation  Defined 
by  Eq.  (  10c) 
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n  Similarity  Coordinate  -  Outer  Layer,  Eq.  (30) 

(x)  Function  Appearing  in  Solution  for  Skin 
Friction  Eq.  (1 15a) 

•j„  Constant  Appearing  iu  Solution  for  Skin 

Friction  Eq.  (115b) 

k;  von  Karmen  Constant  (k  =  0.43) 

<i  Scaled  Von  Karmen  Constant  (ki  :  <a  = 

0.80) 

X  Constant,  \  =  (ae/c;k)  Eq.  (45b) 

A  Local  Skin  Friction  Parameter, 


\t  Turbulent  Viscosity,  Eq.  (1) 

v  Laminar  Viscosity 

p  Density 

ok,a£  Turbulent  Model  Constants  (m.  =  1,  ae  = 
1.3) 

%  Skin  Friction 

Tw  r  Reference  Value  of  Skin  Friction 

Tw  Normalized  Skin  Friction  Tw=iw/TW)r 

Superscripts 

Denotes  Diffusion  Layer  Quantities 
Denotes  Transition  Layer  Quantities 
'  Prime  -  Denotes  Derivative  with  Respect 

to  Similarity  Variables,  n,  C 
*  Denotes  Dimensional  Variable 

+  Wall  Layer  Variable 

(0)  Perturbation  Index 

Subscripts 

D  Diffusion  Layer  Solution 

e  Wall  Set  Edge  Value 

m  Denotes  Value  at  the  Maximum  of  the 

Streamwise  Velocity 
w  Value  at  Wall  (y,  y,  y  =  0) 

o  Perturbation  Index 

cp  "Common  Part” 


1.  Introduction 

In  this  paper  we  examine  the  theoretical  struc 
ture  of  radial  and  planar  wall  jets  spreading  out  on 
a  flat  surface  in  stagnant  surroundings.  Although 
geometrically  simple,  these  flows  raise  many  theo¬ 
retical  questions  relating  to  turbulence  modeling 
and  to  wall  jet  similarity  and  asymptotic  structure. 
Wall  jets  are  of  interest  in  turbulence  modeling  be¬ 
cause  they  combine  in  a  single  flow  field  important 
aspects  of  both  free  and  wall  bounded  turbulence. 
They  are  also  simple  examples  of  a  flow  with  a  ve¬ 
locity  maximum  that  is  not  coincident  with  a  van¬ 
ishing  turbulent  shear  stress;  a  feature  that  can¬ 
not  be  described  by  turbulence  models  employing 
an  eddy  viscosity  hypothesis. 

As  discussed,  for  example,  by  Townsend1 
turbulent  wall  jets  developing  in  stagnant  surround 
ings  belong  to  a  general  class  of  flows  referred  to 
as  self  preserving.  Such  a  flow  is  similar  (or  self- 
preserving)  if  the  form  of  the  solution,  when  non- 
dimensionalized  with  respect  to  certain  velocity 
and  length  scales,  is  independent  of  r,  the  coordi¬ 
nate  parallel  to  the  shear  layer.  These  velocity 
and  length  scales  are  usually  simple  power  law 
functions  of  the  streamwise  coordinate.  The  simi¬ 
larity  transformation  reduces  the  original  partial 
differential  equations  to  ordinary  differential  equa¬ 
tions.  If  the  similarity  is  exact,  the  solutions  of 
the  ordinary  differential  equations  are  also  exact 
solutions  of  the  original  governing  equations.  The 


conditions  that  need  to  be  satisfied  for  similar  so¬ 
lutions  to  exist  in  turbulent  flow  are  discussed  by 
Townsend2,  Rotta-1,  Rodi4  and  many  others.  Exam¬ 
ples  of  similarity  solutions  in  free  turbulent  flows 
are  well  known  (e.g.  jets,  mixing  layers)  and  are 
given  in  standard  fluid  mechanics  texts  (e.g.  Ref. 
5).  Wall  bounded  turbulent  flows  generally  do  not 
satisfy  the  conditions  for  exact  similarity  except  in 
very  special  circumstances  that  are  much  more  re¬ 
strictive  than  in  laminar  flow. 


Wall  bounded  flows,  however,  do  exiubit  a  more 
general  form  of  similarity  referred  to  as  asymptotic 
similarity.  In  this  case,  if  the  external  conditions 
permit  it,  the  boundary  layer  approaches  a  similar 
state  far  downstream.  The  similar  solution  is  then 
the  leading  term  in  an  asymptotic,  coordinate  type 
expansion,  valid  for  r>-  .  In  this  case,  the  leading 
order  similarity  solution  is  generally  not  an  exact 
solution  of  the  original  governing  equations  but  is 
only  an  approximation  that  becomes  more  exact  as 
distance  increases.  Nonsimilar  effects  arise  from 
higher  order  terms  of  the  underlying  asymptotic  ex¬ 
pansions.  As  is  well  known,  turbulent  wall  layers 
develop  in  a  multi  layer  structure  consisting  of  two 
or  more  layers.  The  solution  in  each  layer  can  be 
developed  in  a  far  field  similarity  form ,  but  the  form 
of  the  similarity  solution  and  the  velocity  and  length 
scales  differ  in  each  of  the  layers.  Thus,  the  overall 
solution  never  achieves  a  completely  self- preserving 
asymptotic  state  because  the  solution  in  each  region 
develops  at  a  different  rate.  Nevertheless  the  solu¬ 
tion  in  each  region  correctly  matches  its  neighbors, 
in  the  sense  of  the  method  of  matched  asymptotic  ex¬ 
pansions,  and  the  overall  solution  is  a  rational  asymp¬ 
totic  solution  valid  far  downstream. 

Glauert6  developed  an  approximate  analytic 
solution  for  both  planar  and  radial  wall  jets  in  still 
air.  He  was  able  to  construct  a  solution  in  near¬ 
similarity  form  by  assuming  the  flow  to  be  divided 
into  two  regions:  an  outer  free  jet  type  flow,  and 
an  inner  conventional  turbulent  boundary  layer. 

The  two  solutions  were  matched  at  the  position  of 
the  velocity  maximum  where  the  turbulent  Reynolds 
stresses  were  also  assumed  to  vanish.  The  solution 
in  the  outer  layer  was  obtained  with  a  constant  eddy 
viscosity  model  appropriate  to  free  jet  flows  while 
the  solution  in  the  inner  layer  was  assumed  to  be  the 
same  as  in  turbulent  pipe  flow.  Townsend1  also  de¬ 
veloped  a  near-similar  two  layer  solution  for  the  wall 
jet  but  he  employed  a  standard  logarithmic  wall  layer 
for  the  inner  solution  instead  of  the  power  law  pro¬ 
file  used  by  Glauert  Comparison  of  Glauert's  theory 
with  experiment  shows  that  it  can  accurately  predict 
the  velocity  profiles  and  most  other  features  of  the 
flow.  The  theory  indicates  that  the  inner  and  outer 
layers  grow  at  slightly  different  rates  and  that  the 
exponents  in  the  velocity  and  length  scales  slightly 
differ  from  the  similarity  values  of  free  jets  These 
effects  are  due  to  the  approximate  nature  of  the 
similarity  of  wall  jets. 


Two  controversial  aspects  of  the  theoretical 
solution  arc  the  values  of  the  exponents  in  the  veloci¬ 
ty  and  length  scales  and  the  question  of  the  existence 
of  a  logarithmic  behavior  in  the  velocity  profile  near 
the  wall.  For  a  radial  free  jet,  the  similarity  solu¬ 
tions  indicate  that  the  length  scale,  5j/2>  increases 
linearly  with  radial  distance,  r,  while  the  velocity 
scale,  Urp  decays  as  r'l.  The  experimental  results 
of  Bakke1  and  Porch,  Tsuei  and  Cermak6  for  radial 
wall  jets  are  consistent  with  Glauert’s  prediction  for 
the  velocity  scale,  both  showing  a  10-15%  decrease 
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in  the  exponent  from  the  similarity  value.  However, 
the  same  (lata  also  indicates  a  5  to  10%  decrease  in 
the  exponent  of  the  length  scale  for  which  the 
Glauert  theory  predicts  a  2%  increase  above  the  simi¬ 
larity  value.  For  the  planar  free  jet  the  similarity 
exponent  for  the  velocity  scale  is  -1/2  while  that  for 
the  length  scale  is  1  which  is  the  same  as  the  value 
in  a  radial  jet.  The  data  for  the  velocity  scale  in 
planar  wall  jets  was  summarized  by  Hammond^.  The 
experimental  results  indicated  a  decrease  of  exponent 
in  the  velocity  scale  by  about  5-10%  from  the  similar¬ 
ity  value  which,  as  in  the  radial  case,  is  consistent 
with  Glauert's  theory.  The  results  for  the  length 
scale  in  the  planar  case  are  somewhat  equivocal. 

The  experimental  results  of  Schwarz  and  Cosart*0 
and  Meyers,  Schauer  and  Eustis11  indicate  that  the 
exponent  for  the  length  scale  is  equal  to  the  similar 
ity  value,  1.  while  the  Bradshaw  and  Gee1'-  data 
indicate  a  value  10%  lower. 

There  is  considerable  uncertainty  as  to  whether 
a  conventional  logarithmic  law  of  the  wall  region  de¬ 
velops  in  the  velocity  profile  of  a  wall  jet.  The  same 
theoretical  arguments  that  support  a  log  law  in 
boundary  layers  and  pipe  flows  suggests  that  it 
should  also  arise  in  wall  jets.  However,  the  experi¬ 
mental  results  of  Bradshaw  and  Geel2,  Poreh.  Tsuei 
and  CermakS,  Schwarz  and  CosartlO  and  others,  as 
well  as  the  analysis  of  Hammond9  all  indicate  that  wall 
jets  do  not  exhibit  a  well  defined  log  law  region. 

These  results  indicate  that,  if  a  log  law  profile  exists 
at  all,  it  must  extend  only  over  a  very  narrow  region, 
and  with  constants  that  differ  significantly  from  the 
standard  values  appropriate  to  boundary  layers  and 
pipe  flows.  Nevertheless,  a  conventional  log  law 
behavior  was  assumed  in  setting  the  boundary  con¬ 
ditions  in  recent  calculations  by  l.juboja  and  Rodi  • 3 
using  the  k-<_  model  for  a  wall  jet  on  flat  surface  and 
by  Gibson  and  Younis'4  in  calculations  with  a  Rey 
nolds  stress  model  on  curved  surfaces. 

Wall  jets  are  known  to  spread  less  rapidly  than 
free  jets.  I.juboja  and  Rodi  la  note  that  the  reduced 
spread  rate  of  wall  jets  could  not  be  accounted  for 
with  the  standard  k  .  model  employing  the  same  con 
stants  that  successfully  predict  free  jet  flows.  The 
standard  k model  overpredicts  the  wall  jet  spread 
rate  by  30  percent.  They  suggest  that  this  failure 
is  due  mainly  to  the  inadequacies  of  the  turbulence 
model  with  regard  to  the  "damping"  effect  of  the  wall 
on  the  normal  velocity  fluctuations  rather  than  due 
to  the  direct  action  of  the  wall  shear  stress.  Their 
results,  obtained  with  a  modified  k- c  model  that  ac¬ 
counted  for  the  wall  damping  effect,  showed  much 
improved  predictions  of  the  spread  rate. 

In  order  to  address  the  foregoing  issues  we 
thought  it  desirable  to  develop  a  more  rigorous  ana¬ 
lytical  solution  of  the  wall  jet  problem.  In  the  present 
paper  we  will  not  be  concerned  with  turbulence 
modeling  per  se.  The  turbulence  modeling  aspects 
of  the  wall  jet  problem  have  recently  received  much 
attention  and  the  available  data  and  turbulence 
models  have  been  evaluated  at  the  1980  81  Al-'OSR 
Stanford  conference  on  complex  turbulent  flows  l:’. 
Instead,  in  the  present  paper  we  are  concerned  with 
the  overall  theoretical  structure  of  wall  jet  flows  and 
with  the  development  of  a  general  theoretical  up 
proach  to  this  class  of  flows.  In  order  to  keep  the 
analysis  as  simple  as  possible  and  to  produce  definite 
results  that  can  be  compared  with  available  numeri 
cal  solutions  we’ve  based  the  present  study  on  the 
standard  k  .  model  including  the  wall  damping  modi 
fications  proposed  by  I.juboja  and  Rodi^.  3  he  pres¬ 
ent  theoretical  study  indicates  that  a  key  feature  of 


the  wall  jet  flow  is  that  the  Reynolds  stresses  in  the 
outer  layer  are  an  order  of  magnitude  larger  than  the 
levels  near  the  wall.  The  analysis  to  follow  will  show 
that  this  feature  dominates  the  theoretical  develop¬ 
ment  and  leads  to  a  characteristic  multi-layer  struc¬ 
ture  that  is  typical  for  flows  with  relatively  large 
Reynolds  stress  levels  in  the  outer  region.  We  be¬ 
lieve  that  the  theoretical  approach  developed  here  is 
general  and  may  prove  useful  for  analyzing  other 
flows  with  the  k-s:  or  other  turbulence  models. 

Our  approach  to  the  wall  jet  problem  is  based 
on  a  rational  asymptotic  analysis  of  the  k  _  model 
equations  using  the  method  of  matched  asymptotic 
expansions.  It  is  an  outgrowth  of  the  asymptotic 
methods  developed  for  boundary  layers  and  for  vari 
ous  turbulent  viscid  in  viscid  interaction  problems 
(see  Refs  16.  17  for  a  review  of  this  latter  work). 

In  ail  of  these  works  the  solution  was  developed  as 
an  asymptotic  expansion  in  terms  of  the  Reynolds 
number,  or  equi valcntly  in  terms  of  a  reference 
value  of  the  friction  velocity.  In  these  problems 
the  asymptotic  approach  leads  to  a  multi  layer  struc¬ 
ture  with  the  velocity  in  the  outer  part  of  the  bound¬ 
ary  layer  developing  a  small  defect  form.  The  pres¬ 
ent  work  represents  a  significant  departure  from 
these  studies  in  that  the  velocity  profile  in  the  outer 
region  of  a  wall  jet  is  not  in  a  small  defect  form  Be¬ 
cause  of  this,  the  wall  jet  develops  a  multi-layer 
structure  that  is  considerably  more  complex  than  that 
which  arises  in  a  boundary  layer. 

In  the  present  paper  we  are  primarily  inter¬ 
ested  in  the  solution  in  the  far  field,  similarity 
region.  The  asymptotic  solution  is  carried  out  in 
two  stages.  The  first  is  based  on  a  two  parameter 
expansion  of  the  full  Reynolds  averaged  equations 
with  a  k  ■_  model  of  turbulence  quantities.  One  of 
the  small  parameters,  y,  is  related  to  the  nondimen 
sional  friction  velocity,  u1,  defined  by  the  surface 
shear  stress.  It  is  basically  a  Reynolds  number 
parameter,  y  =  0(ln  Re)"',  that  primarily  controls 
the  shear  stress  effects  induced  by  the  wall.  The 
aforementioned  asymptotic  theories  for  boundary 
layers  and  flow  in  pipes  are  basically  expansions  in 
this  parameter.  The  other  small  parameter  appear 
ing  in  our  analysis.  . .  is  a  measure  of  the  turbu 
lenee  levels  in  the  outer  free  jet  part  of  the  flow. 

For  Reynolds  numbers  approaching  infinity  ,  the  wall 
jet  now  approaches  a  limiting  inviscid  (Reynolds 
number  independent)  turbulent  state  that  is  elosely 
related  to  classical  free  jet  flows.  Since  for  Re  -»•  «> 
there  are  no  other  nondimensional  parameters  appear¬ 
ing  in  the  underlying  deterministic  governing  equa 
tions .  i  must  be  a  basically  statistical  parameter  in 
troduced  by  the  coarse  graining  or  averaging  of  the 
turbulent  flow.  Our  analysis  indicates  that  the 
■  -0  limit  leads  to  a  standard  boundary  layer  de¬ 
scription  with  filling  the  role  ol  the  noudimen 
sional  viscosity  coefficient  (i.e.  lie  ').  The  value 
of  i  is  related  lo  t  u  modeling  constants  arising  in 
the  chosen  turbulence  closure.  n  the  present 
analysis,  in  which  we  employ  the  k  ■  model,  a  is 
identified  with  the  constant,  c  ,  appearing  in  the 
eddy  viscosity  formula  for  the  Reynolds  shear  stress. 
The  constant  c  is  usually  given  by  e  0.051 
which  is  not  impressively  small  for  an  asymptotic 
analysis.  Nevertheless,  the  use  of  a  small  c  as 
ymptotic  approach  is  plausible  since  its  main  effect 
is  to  reduce  the  Reynolds  averaged  equations  to 
boundary  layer  form.  This  is  known  to  lie  a 
generally  useful  approximation  for  both  free  and  wall 
bounded  turbulent  flows,  although  not  nearly  as 
good  as  for  laminar  flows. 
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Thus,  in  tile  first  stage  of  the  analysis  we  ex 
pand  the  k -v  version  of  the  full  Reynolds  averaged 
equations  for  y.  u->-0.  The  a  expansion  reduces  the 
problem  to  a  classical  boundary  layer  formulation  to 
lowest  order.  The  expansion  for  y ->-0  leads  to  a  four 
layer  description  of  the  wall  jet.  The  outer  layer  is 
closely  related  to  a  free  jet  flow  while  the  innermost 
layer  is  a  classical  law  of  the  wall  region.  Two  addi 
tional  intermediate  layers  are  needed  to  effect  the 
matching  of  the  outer  and  inner  layers  and  to  com 
plete  the  solution.  The  analysis  leads  to  a  set  of 
either  partial  or  differential  equations  governing 
the  solution  in  each  region.  The  resulting  descrip 
tion  applies  to  the  complete  development  of  the  wall 
jet  How ,  both  in  the  near  and  far  field ,  except  for 
the  immediate  vicinity  of  an  impingement  region, 
where  the  boundary  layer  approximation  of  constant 
pressure  across  the  shear  layer  is  not  valid. 

In  the  second  stage  of  our  analysis  we  expand 
the  solution  in  each  layer  in  terms  of  distance  along 
the  wall  to  obtain  a  far  field  solution  in  each  region. 
The  leading  order  far  field  solution  in  each  layer  is 
expressible  in  a  similarity  form  which  turns  out  to  be 
an  exact  solution  of  the  equations  governing  the 
solution  in  each  layer. 

One  of  the  interesting  consequences  of  our  ap¬ 
proach  concerns  the  character  of  the  far  field  solu  - 
tion  in  the  outermost  layer  of  the  wall  jet.  Because 
turbulent  shear  stress  levels  in  the  outer  region  are 
large,  0  (a),  compared  to  the  0  (u3)  levels  in  the 
inner  wall  layers,  the  turbulent  shear  stress  in  the 
outer  solution  must  satisfy  a  zero  boundary  condi¬ 
tion  at  the  wall.  This,  together  with  the  fact  that 
the  normal  component  of  velocity  is  also  zero  at  the 
wall  leads  to  an  outer  solution  that  is  similar  in  many 
respects  to  a  symmetrical  free  jet  flow.  Because  the 
turbulent  shear  stress  vanishes  at  the  wall,  the 
momentum  flux  is  conserved  in  the  outer  layer  of  the 
wall  jet.  to  lowest  order.  Consequently  the  outer 
solution  has  the  same  similarity  form  in  the  far  field 
as  a  free  jet  and  the  exponents  in  the  velocity  and 
length  scales  are  equal  to  the  free  jet  similarity 
values.  However,  the  outer  layer  solution  differs 
from  the  free  jet  solution  because  the  presence  of  the 
wall  directly  influences  the  outer  solution  through 
the  wall  boundary  conditions  on  the  turbulent  energy, 
k,  and  dissipation,  _,  leading  to  a  zero  turbulent 
viscosity  at  the  wall.  Because  of  these  boundary 
conditions  the  wal'  jet  profiles  in  the  outer  region 
differ  considerably  from  those  of  a  free  jet  solution, 
particularly  near  the  symmetry  plane.  It  is  also  of 
interest  to  note  that  the  difference  in  the  wall 
boundary  conditions  for  k  and  ■_  considerably  de 
crease  the  spreading  rate  of  the  wall  jet  compared 
to  a  free  jet,  as  is  observed  in  practice. 

The  solution  in  the  innermost  region  is  repre 
sented  by  the  conventional  law  of  the  wall  formula. 
The  overlap  with  its  neighboring  layer  yields  a  con 
ventional  log  law  behavior  with  standard  values  as¬ 
sumed  for  the  associated  constants.  However,  be 
cause  of  the  rapid  variation  of  the  solution  in  the 
other  layers  near  the  wall  the  extent  of  the  log  law 
overlap  is  much  reduced  compared  to  that  in  standard 
turbulent  boundary  layers  and  this  tends  to  support 
the  experimental  findings  reported  previously#,  1  0,1 2. 

In  the  present  work  we've  developed  a  complete 
formulation,  including  the  matching,  for  the  leading 
order  solution  in  each  of  the  four  layers.  We've  re 
dueed  the  governing  equations  in  each  layer  to  far 
field  similarity  form  and  have  obtained  solutions 
for  all  the  resulting  ordinary  differential  equations. 


We  compare  certain  features  of  our  far  field  solu 
tion  with  numerical  solutions  of  the  partial  differen¬ 
tial  equations  for  the  wall  jet  problem  recently  pub 
lished  by  l.juboja  and  Rodi13  using  the  same  form  of 
k  l  turbulence  model.  Wc  also  compare  our  solutions 
for  the  velocity  and  length  scales  with  the  experi¬ 
mental  correlations  for  planar  wall  jets  presented  by 
Narasimha,  Yenga  Naravan  and  Parthasarathy  1 8. 

It  is  of  interest  to  note  that  the  present  theory  leads 
to  the  same  similarity  variables  proposed  by  those 
authors  on  the  basis  of  a  dimensional  analysis  of  the 
data  and  on  the  recognization  that  the  initial  monicn 
turn  flux  is  the  only  relevant  initial  condition  affect¬ 
ing  the  far  field  solution. 

The  solutions  developed  to  date  deal  only  with 
the  leading  order  far  field  similarity  solution  and 
hence  predict  velocity  and  length  scales  that  vary 
in  accordance  with  the  similarity  values  for  the  ex¬ 
ponents.  In  order  to  address  the  nonsimilar  effects 
associated  with  the  shifts  in  velocity  and  length 
scale  exponents  observed  in  the  data  wc  would  need 
to  determine  the  next  terms  in  the  far  field  similarity 
expansions.  Current  efforts  arc  directed  at  the 
computation  of  these  terms. 

II .  Governing  Equations.  Asymptotic  Formulation 
and  structure 

We  assume  the  wall  jet  is  formed,  either  by  im¬ 
pinging  a  jet  onto  a  flat  surface  or  by  blowing- 
through  a  narrow  slot  located  near  the  origin  as 
shown  in  Fig.  1.  The  initial  conditions  of  the  jet  are 
characterized  by  velocity  and  length  scales  uj*  and 
b*,  where  these  are  dimensional  quantities.  All 
other  quantities  referred  to  in  the  paper  are  made 
dimensionless  using  these  scales  and  the  value  of 
density  in  the  initial  jet.  The  quantities  um  and  6m 
are  the  value  and  position  of  the  maximum  velocity 
while  v  i/2  denotes  the  position  of  the  half  velocity 
point.  Our  analysis  is  based  on  the  Reynolds  aver¬ 
aged  equations  using  the  version  of  the  k  -  v  model 
employed  by  Ljuboja  and  Rodi13  in  their  numerical 
study  of  the  planar  wall  jet.  Their  version  includes 
modifications  to  account  for:  the  wall  "damping" 
effect,  and  the  effect  of  the  ratio  of  production  to 
dissipation.  P It,  on  the  turbulent  viscosity  coeffi¬ 
cient,  vt.  We  closely  follow  the  Ljuboja  and  Rodi13 
formulation  and  notation  except  for  a  minor  change 
in  the  continuity  equation  to  allow  for  the  radial  wall 
jet.  In  their  formulation  the  turbulent  viscosity  co¬ 
efficient  is  written  in  the  form 

r?  -  cwk*!4.*  ,,*  (1) 

where  the  function  cp  is  given  by 

cu  =cU0Gu(P/f,  f)  (2) 
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Fig.  1  Wall  Jet  in  Still  Air 
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And  cu0  is  one  of  the  turbulence  model  constants  ^  g<0'U><0 

( taken  to  be  equal  to  c u0  =  0 . 09) ,  G  M  h  G  j  G  2  is  “  "  ' 

defined  in  Ljuboja  and  Rodi  and  f  is  their  wall  "damp-  where  <j  is  a  "reduced"  v 
ing"  function,  defined  by,  <.  :  <  A1/4. 

I  ^*3/2  3 / 4  j  coordinates  deTined  by 

f  =  MIN  <1,  — > 

(  y  e  j  (3)  r*/b*  r  x 


von  Karmen  constant 
.  and,  x,  y  are  stretched 


in  here  k  is  the  von  Karmen  constant . 

The  asymptotic  analysis  is  carried  out  in  terms 
of  two  basic  small  parameters,  a  and  y.  The 
"boundary  layer"  parameter  a  is  defined  to  equal 


while  the  "friction  parameter"  y  is  defined  in  terms 
of  a  reference  value  of  the  nondimensional  friction 
velocity,  uT,  by 

y  =uTir/o<1*i>/'  (5 

where  the  friction  velocity  is  defined  by 

u,  -  V'r„  '  =  \'"Cf72  (6 

From  turbulent  wall  layer  flows  we  know  that  the 
friction  velocity  depends  logarithmically  on  the 
Reynolds  number  Rej, 


where  *  denote  dimensional  quantities.  We  call 
attention  to  the  appearance  of  the  scaling  param¬ 
eter,  a,  in  the  definitions  of  the  dependent  vari¬ 
ables.  The  j  dependence  follows  from  a  considera¬ 
tion  of  the  balance  of  the  overall  momentum  flux . 

We  also  note  that  the  boundary  layer  scaling  appears 
in  the  r  coordinate  direction  with  the  thickness  of 
the  wall  jet  equal  to  the  order  of  magnitude  of  the 
slot  width  or  other  characteristic  dimension  of  the 
impinging  jet.  This  scaling  implies  that  the  solution 
holds  in  a  region  that  is  far  downstream  compared 
to  the  characteristic  dimension  of  the  initiating  jet. 
The  substitution  of  Eqs.  (10)  into  the  full  Reynolds 
averaged  equations,  using  the  above  stretchings, 
and  taking  the  limit  a  •+  0  leads  to  the  standard 
boundary  layer  form  of  the  equations  of  mass,  mo¬ 
mentum,  turbulent  energy  and  dissipation,  written 
in  the  form,  (with  the  perturbation  superscripts 
deleted  for  simplicity) 

au  av  ju 

3x  +  3y+x  (12) 


Re,  -  u  ,*b*/i'*  (  8) 

In  our  approach,  we  seek  solutions  of  the  full 
Reynolds  averaged  equations  in  the  double  limit, 

a  —  0,  y  —  0  (8) 

Typical  numerical  values  for  these  parameters  for 
practical  wall  jet  flows  are  a.  y  ~  0.1.  Since  these 
values  are  not  impressively  small  we  expect  the 
higher  order  terms  to  be  important  in  most  situa¬ 
tions.  Nevertheless,  we  assume  they  are  small 
enough  for  an  asymptotic  approach  to  be  useful. 

The  parameter,  «.  enters  the  formulation  in 
virtually  the  same  manner  as  the  laminar  viscosity 
coefficient,  and  therefore  expansion  of  the  solution 
in  powers  of  i  leads  to  the  boundary  layer  equa¬ 
tions  just  as  in  laminar  flow.  The  a  expansion  takes 
the  form 


U*/u?o"<-U  U,(”<x,  y;  y)  +  0(<v1/J)  +  O  (a)  (10a) 

V7u,V!*‘>/4  ?V  V<0,(x,  y;  y)+0(o‘/!)  +0(o)  .  > 
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The  Reynolds  stress,  g,  is  computed  from  the  mean 
velocity  via  an  eddy  viscosity  formula, 

g=''t|^  .  vt  Guk24  (16 

with  the  eddy  viscosity  given  by  Eq.  (1).  We  use 
standard  values  of  the  model  constants  given  by: 


=  0.09,  ok 


=  1.3  cEj  =  1.44  and  cE2 


’ (x,  y;  y )  +0(o1/s)  +  0(o )  (  10c) 


'(x,  y;  y)  +  0(o1/!)  +  Q(o ) 


1.92.  The  a  expansion  provides  the  formal  basis  for 
using  boundary  layer  theory  for  free  turbulent  flows 
as  well  as  for  the  wall  jet  flow  considered  in  this 
paper.  The  second  term  in  the  a  expansion,  Eq . 
(10),  arises  from  small  pressure  gradients  generated 
in  the  outer  entrainment  flow  field  above  the  wall 
jet.  The  disturbances  in  this  region  arc  induced  by 
the  entrainment  of  external  fluid  into  the  wall  jet. 
The  0  (  i)  terms  are  duo  to  higher  order  terms 
(i.e.  normal  pressure  gradient,  streamwise  diffusion 
etc.)  in  the  Reynolds  averaged  equations  that  are 
not  included  in  the  boundary  layer  equations. 

The  ,  dependence  in  the  solution  indicated  in 
Eq.  (10)  arises  from  the  laminar  shear  stress  term  in 
Eq.  (13),  which  is  important  only  in  the  wall  layer 
near  the  surface.  The  expansion  for  y  >-0  leads  to  a 
solution  with  the  complex  four  layer  structure 
sketched  in  Fig.  2.  consisting  of:  an  outer  free  jet 
type  flow,  a  conventional  equilibrium  wall  layer,  and 
two  additional  layers  that  are  needed  to  effect  a 
matching  between  the  outer  and  inner  layers.  The 
solution  in  the  outer  layer  satisfies  the  exact  bound¬ 
ary  layer  equations.  The  Reynolds  number  term 
appearing  in  the  momentum  equation  is  negligible  to 


all  orders  in  the  outer  expansion  so  that  the  outer 
layer  is  a  fully  turbulent  flow  that  is  not  directly 
influenced  by  the  laminar  viscosity  coefficient.  Of 
course,  there  are  Reynolds  number  effects  in  the 
outer  solution,  but  they  enter  only  through  the  in¬ 
ner  matching  conditions.  The  length  scale  for  the 
outer  region.  Ay  -  0(1),  indicated  in  Fig.  2,  follows 
from  boundary  layer  type  scaling  considerations. 

Since  the  spread  rate  of  the  wall  jet  is  0(  »l/2),  the 
thickness  of  the  outer  layer,  up  to  a  distance  .  r  = 
0(a'l/2)  must  be  the  same  order  as  the  characteris¬ 
tic  dimension  of  the  initial  jet.  Since  the  turbulent 
shear  stress  and  energy  are  an  order  of  magnitude 
larger  than  the  =  0(u?)  levels  at  the  wall, 

the  outer  solution  for  g  and  k  must  vanish  at  the 
wall.  From  the  momentum  equation  it  also  follows 
that  the  outer  solution  for  g  must  have  a  linear  be¬ 
havior  near  the  wall.  Detailed  analysis,  summarized 
in  the  following  section,  indicates  that  the  outer  so¬ 
lutions  for  the  turbulent  energy  and  dissipation 
have  a  non -analytic ,  algebraic  behavior  near  the  wall. 

The  solution  in  the  innermost  layer  is  described 
by  a  conventional  equilibrium  wall  layer  that  is  com¬ 
pletely  determined  by  the  small,  0(  y*^  l+i) value  of 
the  wall  shear  stress.  The  solution  in  the  wall  layer 
has  a  standard  logarithmic  law  of  the  wall  behavior 
at  its  upper  limit,  y+  ;  yuT/v+<«.  Because  of  the 
algebraic  behavior  of  the  outer  solution  for  y->-0.  the 
outer  and  wall  layer  solutions  do  not  match  and  ad¬ 
ditional  intermediate  layers  are  required.  Two  ad¬ 
ditional  layers  are  needed:  a  diffusion  layer,  and  a 
transition  layer.  The  thickness  of  the  diffusion  lay¬ 
er  is  determined  by  the  condition  that  the  linearly 
varying  part  of  the  solution  for  the  shear  stress  is 
of  the  same  order  of  magnitude  as  the  value  of  the 
wall  shear  stress.  This  leads  to  a  thickness  y  = 

0( y 2)  as  indicated  in  the  figure.  This  is  the  region 


in  which  the  velocity  has  a  maximum  and  the  turbu¬ 
lent  shear  stress  has  a  zero.  We  call  it  a  diffusion 
layer  because  of  the  dominant  influence  of  turbulent 
diffusion  in  this  layer.  Detailed  analysis  indicates 
that  the  dissipation  is  balanced  entirely  by  turbulent 
diffusion  in  this  region.  The  dominant  terms  in  the 
turbulent  energy  and  diffusion  equations  in  each 
layer  are  indicated  in  Fig.  2  by  the  symbols: 

A  (advection),  P  (production),  Dp  (diffusion),  and 
D .  (dissipation) . 

The  solution  of  the  diffusion  layer  equations 
also  exhibits  an  algebraic  behavior  at  the  wall  which 
does  not  match  the  logarithmic  behavior  of  the  wall 
layer  solution.  This  requires  the  introduction  of 
the  transition  layer  indicated  in  Fig.  2.  The  length 
scale.  Ay  =  y^a,  for  this  layer  is  determined  by 
the  requirement  that  it  matches  to  both  the  diffusion 
and  wall  layer  solutions.  The  constant  "a"  is  a 
function  of  the  turbulence  modeling  constants, 
which  for  the  standard  values  listed  previously  is 
given  by .  a  as  0.  27, 

Ill.  Outer  Region 

The  solution  in  the  outer  region  is  determined 
by  expanding  the  boundary  layer  equations,  Eqs 
12-16,  in  terms  of  the  friction  parameter,  Y  .  Since 
the  viscous  term  appearing  in  the  momentum  equa¬ 
tion  is  exponentially  small  in  Y  it  can  be  neglected 
to  all  orders.  Thus,  the  leading  order  term  in  the 
outer  solution  satisfies  the  exact  boundary  layer 
equations  with  the  laminar  shear  stress  term  set  to 
zero.  We  assume  that  the  solution  in  the  outer  re¬ 
gion  can  be  expanded  in  the  form . 

U(x,  y;  >)  -  u0(x,  y;  <))  +  •••  (17) 
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Fig.  2  Wall  Jet  Flow  Field  Structure 
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;md ,  for  simplicity,  we  drop  the  perturbation  sub 
scripts  in  the  subsequent  developments. 


For  a  stagnant  outer  stream  the  k  -  model 
leads  to  a  shear  flow  with  a  definite  outer  edge, 
y  =  ye-  (e.g.  see  Ref.  19)  where  the  solution  sat¬ 
isfies  homogeneous  boundary  conditions 


k  :  c  -  I)  AX  y 


(22) 


The  boundary  conditions  on  the  inner  edge  of  the 
outer  solution  are  determined  by  matching  to  the 
inner  layers.  Since  g  and  k  are  0  (X  2)  at  the  wall 
(and  are  also  small  in  the  other  inner  layers)  the 
outer  solution  for  these  quantities  must  satisfy  zero 
boundary  conditions  at  the  wall.  The  boundary  con 
dition  on  the  dissipaton ,  is  more  difficult  to  fix 
since  it  cannot  be  determined  from  simple  scaling 
arguments  if  it  should  be  equal  to  zero,  infinite  or 
a  finite  value  at  the  wall. 


A  detailed  analysis  of  the  limiting  behavior  of 
the  outer  solution  for  y  -0  indicates  that  matching 
is  possible  only  if  the  dissipaton  in  the  outer  solu 
tion  is  unbounded  at  the  wall.  The  final  inner 
boundary  condition  is  fixed  by  the  requirement  that 
the  normal  component  of  velocity  bo  zero  at  the  wall. 
Thus  the  inner  boundary  conditions  are 


v  ■  g  =  k  -  0 

AT  y  -  0 

( 23a) 

f  —  or 

AT  y  =  0 

(23b) 

t  he  boundary  conditions  given  in  Kqs.  (22)  and 
(23)  are  sufficient  to  completely  determine  the  so¬ 
lution.  As  in  the  free  jet  problem,  an  additional 
wall  boundary  condition  on  u  is  not  required  be¬ 
cause  the  wall  value  of  u  is  determined  from  an  over¬ 
all  momentum  balance.  Since  the  wall  shear  and  nor¬ 
mal  component  of  velocity  are  zero  at  the  wall  the 
momentum  flux  integral,  J.  is  constant  to  lowest 
order 


J  (2 


2-x)J  j  'u!(x,  y;  0)  dy  +  CXy2)  (  24) 


where  the  0  ( y *)  term  arises  from  the  nonzero  value 
of  the  wall  shear  stress.  The  value  of  the  momentum 
flux  is  set  by  the  initial  conditions  in  the  impinge¬ 
ment  region. 


The  above  boundary  conditions  completely  define  the 
outer  solution,  independent  of  the  solution  in  the 
inner  regions.  However,  the  partial  differential 
equations  governing  the  outer  solution  must  be  inte¬ 
grated  numerically  as  an  initial  value  problem  for 
each  prescribed  value  of  the  momentum  flux,  J.  and 
set  of  initial  profiles.  The  solution  in  the  region, 

\r*  =  0  (b*n  ^"1  under  consideration  depends  on 
the  details  of  the  initial  profile  shapes,  which  can 
only  be  determined  from  considcratio-’  of  the  detailed 
solution  in  the  impingement  region. 


Far  Field  Similarity 


constant  to  lowest  order  in  y .  and  consequently ,  the 
outer  solution  behaves  very  much  like  a  turbulent 
free  jet.  In  particular,  the  outer  solution  for  the 
wall  jet  admits  a  similarity  structure  that  is  of  the 
same  form  as  that  for  a  free  jet.  Thus  the  lowest 
order  outer  solution  for  the  streamwisc  and  normal 
components  of  velocity,  turbulent  energy  and  dissi¬ 
pation  can  be  represented  in  the  form, 


uw(x)F'(rj)  <F'«>) 

-  1) 

(25) 

ujx)  (tj  F'(tj)  -  2i_1 

F(>!)) 

(26) 

g  Ui(x)G(b) 

(27) 

k  ujj(x)  K  (v) 

(28) 

(  u^(x)  x'1  H(v) 

(29) 

where  the  primes  denote  derivatives  with  respect  to 
the  similarity  variable ,  r, ,  defined  by 


Because  the  leading  order  solution  for  the  shear 
stress  vanishes  at  the  wall,  the  momentum  flux  is 


y'x 


(30) 


and,  as  defined,  uw(x)  is  the  value  of  u  at  the  wall. 
In  all  that  follows  the  subscript  "w"  denotes  the  wall 
value  (y  =  0,  of  the  associated  variable.  The  func¬ 
tion  uw(x)  xi+1  F(')  is  the  streamfunction  of  the 
leading  order  outer  solution.  The  wall  velocity, 
uw(x),  can  be  determined  from  the  momentum  flux 
constant,  J,  by  rearranging  Eq .  (24)  into  the 
similarity  form 


u„(x)  -  Cj/3  x'(i*‘ "  2 
where  Cj  is  a  constant  defined  by, 

J 


(31) 


1? 


(277)^ 

F2(T))  drj 


(32a) 

(32b) 


The  similarity  equations  for  momentum,  turbulent 
energy  and  dissipation  are 


2i*‘(F2  +  F  F")  F"  ^  =0  (33) 

2]F'K  +  2~FK^K'V  F"2)-  H  =  0  (34) 


|(5  +  3j)F'H  *  2)-*FH' 


ct1(KGu  F 


2) 


(35) 


The  continuity  equation  is  automatically  satisfied 
through  the  introduction  of  the  "streamfunction", 
F  (n).  The  momentum  equation,  Eq .  (33)  can  be 
integrated  once,  yielding 


2>-‘FF'  +  ^12^  f"  -  a 


(36) 


where  the  wall  boundary  conditions  require  the  con¬ 
stant  of  integration  to  vanish.  The  eddy  viscosity 
formula  for  the  shear  stress  takes  the  form. 


GU(K2/H)  F" 


(37) 


Results  from  the  numerical  solutions  of  the  similarity 
equations  will  be  presented  later  in  the  paper.  Al¬ 
though  the  similarity  solution  is  an  exact  solution  of 
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the  full  nonlinear  boundary  layer  equations,  it  is  not 
an  exact  similarity  as  defined  in  the  introduction, 
because  the  wall  matching  conditions  require  the 
introduction  of  additional  higher  order  terms  in  the 
outer  expansion.  Thus  the  similarity  solution  des¬ 
cribed  above  represents  only  the  leading  term  of  a 
far  field  expansion  valid  for  x->  » ,  which  is  an  ex¬ 
ample  of  what  wc  have  called  an  asymptotic  similarity 
in  the  introduction. 

1. hinting  behavior 

The  asymptotic  flow  field  structure  is  strongly 
influenced  by  the  limit  behavior  of  the  outer  solution 
for  y-»  0.  A  conventional  equilibrium  wall  layer  is 
present  near  the  wall  so  that  the  outer  solutions, 
however  many  layers  there  may  be,  must  eventually 
match  to  a  conventional  log  law  for  y-+0.  The  match¬ 
ing  will  lead  to  a  logarithmic  skin  friction  law  of  the 
same  type  ns  arises  in  boundary  layers  and  pipe 
flows.  Since  the  turbulent  shear  stress  and  energy 
are  small  in  the  wall  layer  then  the  outer  solutions 
for  g  and  k  must  vanish  at  the  wall  as  required  by 
the  boundary  condition,  Eq.  (23). 

The  determination  of  the  correct  boundary  con¬ 
dition  to  impose  on  the  dissipation,  e.  proved  con¬ 
siderably  more  difficult  because  it  was  not  possible 
to  deduce  the  proper  condition  thru  a  simple  scaling 
analysis.  The  logarithmic  law  of  the  wall  requires 
the  dissipation  to  behave  like. 


Thus,  the  magnitude  of  .  near  the  wall  depends  on 
the  length  scale  of  the  inner  layers  and  could  be: 
small,  order  one.  or  unbounded  for  y-0.  depending 
on  the  thickness  and  structural  form  of  the  inner 
layers.  However,  the  structure  of  the  flow  near 
the  wall,  in  turn,  depends  critically  on  the  boundary 
condition  imposed  on  and  cannot  be  determined 
without  a  knowledge  of  this  condition.  Thus  the 
correct  wall  boundary  conditon  on  c  and  the  structure 
of  the  solution  near  the  wall  can  only  be  determined 
thru  a  detailed  analysis  of  the  limiting  behavior  of 
the  outer  solution  for  y->-0  for  all  possible  boundary 
conditions  on  We  first  examined  the  consequences 
of  setting  z  to  zero  at  the  wall.  By  trial  and  error 
we  found  that  the  only  solution  satisfying  this 
condition  has  the  form. 


■c„tx)y,  k  —  Atxly,  (  —  Btx)y1/:  <;{<)) 


where  cw(x)  is 


and  A(x)  an  B(x)  are  functions  determined  from  a 
balance  of  terms  in  the  turbulent  energy'  and  dissi¬ 
pation  equations  To  lowest  order,  the  dominant 
terms  are  provided  by  a  balance  of  turbulent  diffu¬ 
sion  ,  production  and  dissipation  which  loads  to. 


uw  v”;~rcJ  w 
vV 

L  ,Tt  (ce:  -  Ce  )  J 


where  (i  w  is  the  limiting  value  of  the  turbulent  vis¬ 
cosity  function  at  the  wall.  Since  the  wall  damping 
function,  f.  approaches  one  and  the  production  - 
dissipation  ratio.  Vi  .  approaches  a  definite  limit  for 
y-»0.  G  w  is  given  by 

G  |i W  'y-j  t  P  (  )w,  1  !  (13) 


The  functions  A(x)  and  B(x)  depend  on  the  model 
constants,  for  the  standard  values  given  previ¬ 
ously  ,  as  well  as  other  sensible  choices  we  find  that 
A(x)  is  strictly  positive  but  that 


which  demonstrates  that  solutions  of  this  type,  with 
■_  equal  to  zero  at  the  wall,  are  not  possible. 

We  next  considered  the  possibility  that  the  dis¬ 
sipation  was  either  unbounded  or  equal  to  a  nonzero 
constant  at  the  wall.  We  assumed  that  the  velocity 
and  turbulent  shear  stress  had  the  same  form  as  in 
Eq  (39)  but  that  the  turbulent  energy  and  dissipa¬ 
tion  had  a  more  general  power  law  variation  lor 
y-*0.  Thus  we  assumed  the  outer  solution  had  the 
following  behavior 

u  g  — c„(x)y,  k  -  Atxly*, 

t  -  B(x)  y'b  FOB  v  -  (I 

where  a  >  0  in  order  to  satisfy  the  boundary  condi¬ 
tion  on  k  (k-»0)  and  b  >  0  in  order  to  avoid  the  con¬ 
sequences  of  t—0  discussed  above.  A  solution  of 
this  type  proved  possible  with  the  dominant  turbu¬ 
lence  terms  involving  a  balance  of  diffusion  and 
dissipation.  Turbulent  production  proved  negligible 
compared  to  both  these  terms  in  this  type  solution. 
The  balance  of  the  diffusion  and  dissipation  terms 
resulted  in  the  following  expressions  for  the  ex¬ 
ponents,  a  and  b. 


b  1  -  (3  2)  a 

-  7  -  (1  -  24X)1/: 
3  5(X  -  2) 


.  X  <"«  ,Tk>  <-■(;  (45b) 


where  the  positive  root  in  Eq .  (45b)  was  chosen  to 
satisfy  the  wall  boundary  condition  on  k.  For  the 
values  of  model  constants  used  here  a  -  0.270  and 
b  =  0.  595.  Since  b  is  positive  the  outer  solution  for 
the  dissipation  is  unbounded  at  the  wall.  A  solution 
with  a  bounded  wall  value  of  .  proved  to  be  impos¬ 
sible.  The  functions  A  (x),  B  (x)  are  to  be  extract 
ed  from  the  numerical  solution  of  the  outer  eq ua 
tions  for  y-*-0.  The  analysis  for  y->-0  shows  that  they 
are  related  by. 

A3/:  B  -  a"1  (2ck  3GU„)1,;  (40) 

The  next  term  in  the  expansion  for  u  near  the  wall 
is  determined  from  an  integration  of  the  eddy  vis¬ 
cosity  equation  to  yield. 


u  —  u„i\)  - 


i  tA"  G  )(2  -  a) 


where  g  0  is  the  limiting  value  of  the  viscosity 
function  at  the  wall  defined  by. 

G„(,  -GM{(P'e)  =  (I,  lj-  (47b) 

The  boundary  conditions  Eqs.  (22)  and  (23) 
completely  fix  the  outer  solution,  independent  of  the 
solution  in  the  inner  layers.  The  general  character 
of  the  outer  solution  for  the  streamwise  velocity,  u. 
is  as  sketched  in  Fig.  2.  The  velocity  increases 
monotonically  from  zero  at  the  upper  edge  of  the  jet 
to  a  maximum.  uw  at  the  wall.  At  the  upper  edge 
the  solution  has  an  algebraic  singularity  of  the  same 
lorm  as  in  a  free  iet'9,  ,.\t  tlx*  wall  the  solution  has 
the  singular  behavior  indicated  in  Eqs.  (II)  and 
(  47).  Thus  the  outer  solution  clearly  cannot  match 
the  logarithmic  behavior  ol  the  wall  layer,  and  as  a 
consequence,  the  skill  friction  cannot  be  determined. 
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IV.  Diffusion  Layer 


The  analysis  of  the  previous  section  indicated 
that  the  outer  solution  exhibits  an  algebraic  behav¬ 
ior  at  the  wall,  in  which  the  turbulent  shear  stress 
and  energy  vanished  while  the  dissipation  becomes 
unbounded  and  the  velocity  approaches  a  definite 
limit,  uw(x).  Because  of  the  algebraic  behavior  for 
y  '0  it  is  not  possible  to  match  the  outer  solution  to 
the  logarithmic  law  of  the  wall.  This  implies  that 
additional  layers  will  be  required  to  achieve  a  match¬ 
ing  to  the  wall  layer  solution.  The  limiting  behavior 
of  the  outer  solution  indicates  that  production  and 
advection  of  turbulent  energy  arc  small  for  y  +0  re¬ 
sulting  in  a  balance  of  diffusion  and  dissipation  at 
the  wall.  This  suggests  that  we  attempt  to  develop 
a  solution  in  a  thin  inner  region  in  which  the  tur 
bulence  characteristics  are  dominanted  by  a  balance 
of  diffusion  and  dissipation.  Because  of  the  impor¬ 
tance  of  turbulent  diffusion  in  this  region  we  call 
it  the  diffusion  layer. 


The  structure  of  the  solution  in  the  diffusion 
layer  is  completely  determined  by  the  behavior  of 
the  turbulent  shear  stress  near  t lie  wall.  In  the  out 
er  solution  the  shear  stress  is  negative  with  a  linear 
decay  to  zero  at  the  wall.  However  the  turbulent 
shear  stress  is  a  positive  constant  in  the  wall  layer 
equal  to  the  small,  0  (Y2u(1+j)/2)  value  of  the  skin 
friction,  tw.  The  outer  solution  breaks  down  where 
the  linearly  varying  shear  stress  becomes  the  same 
order  as  the  wall  shear  stress.  An  analysis,  not  pre¬ 
sented  here,  shows  that  the  breakdown  manifests 
itself  through  the  appearance  of  a  singular  contribu¬ 
tion  in  the  second-order  outer  solution  for  the 
streamwise  velocity.  The  condition  that  the  linearly 
varying  outer  solution  for  the  shear  stress  is  the 
same  order  as  the  wall  shear  estabishes  the  thick¬ 
ness  scale  of  the  diffusion  layer  to  be  A  y  =  n  (7  2). 
Thus,  the  solution  in  the  diffusion  layer  is  express¬ 
ed  in  terms  of  the  stretched  variable,  y,  defined  by 


y  =  r2y 


C  48) 


Appropriate  scalings  for  all  dependent  variables 
follow  from  the  substitution  of  the  stretched  variable 
defined  in  Eq.  (48)  into  the  inner  limit  of  the  outer 
solution  given  by  Eqs.  (44)  -  (47).  Thus  we  are 
led  to  solutions  of  the  form 


U  =  u„  +  y2‘au(x,  y) 


(49) 


g  T2(cwy  +  r„) 


Tw.r[cwy  +  TJ  + 


(50) 


k  ra[Mx,y)|  +  "‘I 


(51) 


y"2b[e(x,  y)  +  •••! 


(52) 


where  cw(x)  is  the  function  of  uw  defined  in  Eq. 
(40)  and  %  is  equal  to  the  ratio  of  wall  shear  to 
its  reference  value. 


-y,/  I  y 


(53) 


Equation  (50)  follows  from  an  integration  of  the  mo¬ 
mentum  equation,  written  in  stretched  variables. 
With  the  scalings  given  above  the  turbulent  energy 
and  dissipation  equations  reduce,  to  lowest  order, 
to  a  balance  of  turbulent  diffusion  and  dissipation; 
a  balance  that  is  consistent  with  the  limiting  be¬ 
havior  of  the  outer  solution.  These  equations  re¬ 
duce  to  second-order  ordinary  differential  equations 


■  ’  '  ^  -1 


for  k  arid  7.  The  only  solutions  that  satisfy  these 
equations  and  match  the  outer  solution  are  given  by. 


A(x)  ya 
B(x)y‘b 


(54) 


(55) 


which  are  precisely  the  limiting  forms  of  the  outer 
solution  with  the  functions  A(x),  B(x)  determined 
from  the  outer  solution  and  the  constants  a,b  given 
by  Eqs  (45).  The  solution  for  the  streamwise  velo¬ 
city  is  found  by  substituting  the  above  expressions 
for  g.  k  and  e  into  the  eddy  viscosity  formula  which 
yields, 

a“  B  .  -  - 

(56) 


dy  A  Guo 


(cwy  *  rjy- 


where  Gy0  is  the  zero  production  limit  of  the  tur¬ 
bulent  viscosity  function  defined  in  Eq  (47b) 


The  integration  of  Eq  (56)  to  obtain  u  introduces  an 
arbitrary  function  of  x  into  the  solution  that  can  be 
evaluated  by  matching  to  the  outer  solution.  The 
matching  requires  this  function  to  vanish  so  that 
the  solution  for  u  can  be  written  in  the  form. 


U/Uy 


[("D 


<y/y  m)'l/2  + ::  <y/ym)''l/:i 


]  [1  -  a  I"1 


(57) 


where  the  functions  um  and  ym  are, 


u  =( _ \  *u-a/2>  (58a) 

m  \a(2  -  a)  A2  C.  /  ym 


“0 


y  m  7w/cw 


(58b) 


since  B  and  tw  are  positive,  cw  is  negative  and 
a  « 0.27,  the  function  um  is  negative  and  yn,  is 
postive.  It  follows  that  the  corresponding  two  term 
expansion  for  the  streamwise  velocity  has  a  maximum 
located  in  the  diffusion  layer  at  y  =  Y  2  ym  where 


the  maximum  value,  Um>  is 


Um  -  uw  +  y2 


(59) 


The  second  term  in  Eq.  (57)  can  be  shown  to  match 
the  outer  solution  given  in  Eq .  (47).  The  overall 
behavior  of  the  solution  for  the  velocity  profile  in 
the  diffusion  layer  is  as  sketched  in  Fig.  2. 


Far  Field  Similarity 


The  diffusion  layer  solution  can  be  placed  in 
similarity  form  by  taking  note  of  the  similarity  forms 
of  the  function  A(x),  B(x),  um(x),  cw(x)  that 

follow  from  the  definition  of  cw  in  Eq.  (40)  and  sim¬ 
ilarity  solutions  given  in  Eqs  (25)  -  (31).  Thus  we 
find 


A(x)  A,u;(x)  x‘a  AsCj  x',,*)*a’ 


(60a) 


B(x)  B,uJ,(x)xb 


B,t'J/2.x'(3/2,(,*l*a'  (60b) 


cw(x)  -il  2X1  •  j)  a^(x)  x" 
-  (1  21(1  •  j)  i1, 


(60c) 


where  As.  Ms  are  constants  determined  from  the 


outer  similarity  solution  and  uw(.<)  is  the  similarity 


•,y-> 


J'J 


-‘  V  ‘-x 


i 


-4 


,V.' 


£ 


form  given  in  Eq.  (31).  If  we  define  the  similarity 
variable  for  the  diffusion  layer  by 


h  =  y/yB 


(61) 


the  solution  for  the  velocity  can  be  written  in  the 
similar  form, 

u/um  -  [(i  -  ^fa/2  -  §  T)‘‘a/2]  [l  -  al'1  <62) 


where,  now, 


ym/x  e  A,(x)  =  2(j  4  l)’1  (tw/i£)  (63a) 


“m/\ 


2(1  +  j)l/2B,  1 

2l*1 

L  a(2  -  a)A2GHJ 

L^J 

l-a/2 


(63b) 


and  where  the  term  in  square  brackets  is  a  constant. 
The  x  dependence  of  the  functions  Al  (x)  and 
A 2  (x)  appearing  in  Eqs  (63)  depends  on  the  solu¬ 
tion  for  the  wall  shear  stress  which  will  be  obtained 
in  the  next  section.  The  far  field  solution  for  the 
maximum  velocity  is  given  in  terms  of  A  2  (x)  by. 


Ura  =  m,(x)  [1  +  y2'aA2(x)] 


(64) 


The  solutions  for  the  turbulent  energy  and  dissipa¬ 
tion  can  be  arranged  into  the  following  similarity 
form , 


k/u2,  =  A3A1(x)  ija 

eymA£  =  BsA,(x)‘-S‘b 


(65) 

(66) 


Limiting  Behavior 


The  diffusion  layer  solution  for  k  and  s  is 
given  entirely  by  the  terms  appearing  in  the  limit¬ 
ing  behavior  of  the  outer  solution  for  y-+0,  so  that 
the  matching  of  these  quantities  is  obvious.  To 
establish  a  match  between  the  outer  and  diffusion 
layer  solutions  for  the  velocity  we  rewrite  the  dif¬ 


fusion  layer  solution, 
obtain 


UD ,  in  outer  variables  to 


Un 


--  -  ’■’(-Jt:) 


FOR  y  -0 


(67) 


where  the  subscripts  "D"  and  "cp”  denote  the  dif¬ 
fusion  layer  solution  and  the  "common  part",  the 
latter  being  given  by  the  inner  limit  of  the  outer 
solution,  Eq  (47).  The  first  term  in  Eq.  (67) 
matches  a  corresponding  term  in  the  outer  solution 
while  the  second  term  represents  an  unbalanced  con¬ 
tribution  that  requires  the  introduction  of  an  0  (y2) 
perturbation  to  the  outer  solution. 

The  diffusion  layer  exhibits  an  algebraic  be¬ 
havior  at  the  wall  that  is  of  the  same  form  as  the 
limiting  behavior  of  the  outer  solution,  except  for 
the  velocity  which  has  the  general  form 


UD  =  ig,(x) 


4  )  y -»/2 

\aA2Gy  /  y 


FORy-O  -68a) 


or,  for  x->  «>  ,  the  similarity  form  given  by 


o.(x)  <1  -  y  '• 


^(idf£-)(x)ai 


FOR  y  -  0 


(68b) 


Because  of  the  algebraic  singularity  in  Eqs  (68) 
the  diffusion  layer  obviously  cannot  mutch  the  wall 
layer  solution.  Thus  an  additional  layer  will  be  re¬ 
quired  to  effect  a  transition  between  the  algebraic 
behavior  of  the  diffusion  layer  solution  and  the 
logarithmic  behavior  or  the  wall  layer  solution. 

V.  Transition  Layer  and  Skin  Friction 

The  thickness  of  the  transition  layer  is  deter¬ 
mined  from  a  consideration  of  the  behavior  the  tur¬ 
bulent  energy  near  the  wall.  The  solution  for  the 
turbulent  energy  in  the  diffusion  layer  indicates 
that  it  approaches  zero  like  k — *-Aya  for  y — ►O. 
Since  k  is  0  (y2)  in  the  wall  layer,  the  diffusion 
and  wall  layer  solutions  for  the  turbulent  energy 
will  be  of  the  same  order  of  magnitude  in  a  layer 
with  a  thickness,  Ay  =  0 (y2/a).  Since  a  ^  0.27, 
the  transition  layer  is  considerably  thinner  than 
the  diffusion  layer  and  is  thicker  than  the  expo¬ 
nentially  thin  (in  terms  of  y)  wall  layer.  The  so¬ 
lution  in  the  transition  layer  is  expressed  in  the 
stretched  variable,  "y,  defined  by 


y  -  y ^  /  a  y 


(69) 


The  appropriate  scalings  for  the  dependent  vari¬ 
ables  in  the  transition  layer  follow  from  the  substi 
tution  of  the  stretched  variable  defined  in  Eq.  (69) 
into  the  inner  limit  of  the  diffusion  layer  solution. 
Thus  we  seek  solutions  in  the  transition  layer  in  the 
form . 


g  =  72  P»(x) 


0(y(2/a,’2)l 

k  =  y2[k(x,  y)  +  •  •  -1 

e  =  y-2b/a[i(x,y>  +  •••  I 
U  =  u„(x)  +  y  u(x,  y)  +  •  •  • 


(70) 

(71) 

(72) 

(73) 


Substitution  of  these  expansions  into  the  turbulent 
energy .  dissipation  and  eddy  viscosity  equations 
yields  the  following  system  of  ordinary  differential 
equations , 


(k2  -r„c;1)  (f~/k2) 


9  /g„  k^_  9k  \ 

9y\(7k  t  Sy ) 

Hi 4®)-^-^-°^’® 


9u 

0y 


k2Gu 


(74) 


(75) 


(76) 


Equation  (70)  indicates  that  the  turbulent  shear 
stress  is  constant  across  the  transition  layer  to  low¬ 
est  order.  Then,  by  matching  to  the  wall  layer  so¬ 
lution  it  is  established  that  the  lowest  order  term  for 
g  is  equal  to  the  skin  friction,  Tw,  as  nornmlized  in 
Eq.  (53).  The  wall  damping  function,  f,  is  also 
equal  to  one  across  the  transition  region,  as  in  the 
diffusion  layer,  but  here  the  production  dissipation 
ratio,  P/t_  ,  is  variable  to  lowest  order.  Thus,  the 
viscosity  function.  Gu,  is  determined  from  the  im¬ 
plicit  relation, 


Gu  =  Gu{P/e,  1] 
P/e  -  G2/k2)G:' 


(77) 

(78) 


Equations  (74)  (78)  are  to  be  solved  subject 

to  boundary  conditions  determined  by  matching  the 
solution  to  the  diffusion  layer  for  y — ►  «.  and  to  the 
law  of  the  wall  for  v — >0.  That  the  transition  layer 
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equations  admit  a  log  law  behavior  for  y — ►()  follows 
from  the  fact  that  all  three  turbulence  source  terms: 
production,  dissipation,  and  diffusion  appear  in  the 
lowest  order  description.  It  is  the  balance  of  these 
terms  in  the  k  -  <_  model  that  leads  to  the  standard 
log  law  at  the  wall.  In  order  to  the  match  to  the 
diffusion  layer  the  solution  must  also  satisfy  the 
following  asymptotic  conditions. 


G„  ^  GuwG(P/e,  1) 
P/e  =  (K2G)"‘ 


(86a) 

(86b) 


where  GyW  is  the  value  of  Gy  at  the  wall,  so  that 
G  is  equal  to  one  for  y=  0.  In  terms  of  this  trans¬ 
formation  the  turbulent  energy  and  dissipation 
equations  are, 


k  -Ay1 

FOR  y  —  <*• 

(79a) 

K"  -  GK2  =  -  1 

(87) 

e  -By"6 

FOR  y  —  «■ 

(79b) 

E"  -  p(K)  E  -  0 

(88a) 

u  (2B~T*  ' 

)  y 'a/2  for  y  —  x 

where  2  (K)  is  defined  by 

VaA2Gu0/ 

(79c) 

l3(K)  -  (rre/<rk)  (c62  GK  -  c(1K"') 

(88b) 

In  order  to  complete  the  solution  and  to  deter 
mine  the  skin  friction  the  transition  layer  must 
match  to  the  logarithmic  law  of  the  wall  for  y — ►O. 
In  the  present  notation  this  leads  to  the  following 
matching  conditions, 


u  —  ( \rj^  /*,)  (lnyc(x))  FOR  y  —  0 


(80a) 


G 


■1/2 


e  —  y 


(80b) 

(80c) 


where  Tw  is  the  local  value  of  the  normalized  skin 
friction,  c(x)  is  n  function  to  be  extracted  from  the 
numerical  solution  of  the  transition  layer  equations 
and  c]  is  a  constant  related  to  the  standard  von 
Karmen  constant,  =  0.43,  by  the  relation 


(81) 


which  tor 


0.09  results  in  a  value  kj  =  0.79 


The  six  conditions  given  in  Eqs.  (79)  and  (80)  are 
sufficient  to  completely  solve  the  fifth  order  system 
of  ordinary  differential  equations  given  by  Eqs. 

(74)  -  (76)  and  to  determine  the  solution  for  the 
skin  friction.  However  the  system  of  equations  is 
highly  nonlinear  and  must  be  solved  by  numerical 
methods. 


Similarity. 


Equations  (74)  -  (76)  can  be  considerably 
simplified  by  first  switching  independent  variables 
from  y  to  u  and  then  placing  the  resulting  equa¬ 
tions  in  similarity  form  through  the  following  scale 
transformation . 


.'o 

'-I 

ii 

7 r 

$ 

(82) 

k  -  kw  K(?) 

(83) 

lyl 

€  -  G  E(?) 

(84) 

Y(?) 


(85) 


wnere  kw  is  the  value  of  k  at  the  wall,  as  given  by 
the  law  of  wall  matching  condition,  Eq.  (80b)  and  e 
is  an  arbitrary  constant,  to  be  determined  by 
matching  to  the  diffusion  layer  solution  for  y  *■  »  . 
In  the  following  will  be  taken  us  the  independent 
variable  which  is  strictly  positive  and  ranges  from 
»at  the  wall  (y=  0)  and  to  r  =  0  in  the  far  field 
(y-+  »).  The  turbulent  viscosity  function  is  ex¬ 
pressed  in  the  form, 


The  transformation  for  Y  (  c)  is  obtained  from  an 
integration  of  Eq .  ( 76) , 


Y  =  J  (GK2/E)d?  (89) 

where  we  have  imposed  the  boundary  condition , 
y  ->■  oo  ( u  *>  “  “>  )  for  Y  -*  0 
in  order  to  match  the  law  of  the  wall. 


The  transformation  removes  all  parameters 
from  the  problem  definition,  aside  from  turbulence 
model  constants,  so  that  K  ( and  E  ( r)  are  univer¬ 
sal  functions  that  need  be  computed  just  once  for 
each  set  of  chosen  model  constants.  More  impor¬ 
tantly  ,  the  transformation  completely  decouples  the 
nonlinear  equation  for  the  turbulent  energy,  K  ( O  , 
from  the  second  order  linear  equation  governing 
E  (r,).  Thus  the  solution  can  be  obtained  sequen¬ 
tially  by  first  solving  Eq  (87)  for  K  (c)  ,  next  Eq 
(88)  for  E  (r)  and  finally  Eq  (89)  for  Y  (  c) .  The 
boundary  conditions  on  K  ( 0  are  determined  by 
matching  to  the  diffusion  layer  solution  for  c,  -*•  0 
(y  -*■  »)  and  to  the  wall  layer  solution  for  Y  oo 
(y  0) .  Detailed  analysis  indicates  that  the  far 
field  matching  condition  can  be  satisfied  by, 


K  (?)  —  <n 


FOR  £  —  0  (y  —  °c) 


(90) 


Because  of  the  chosen  scaling  on  k,  the  boundary 
condition  at  the  wall  is  given  by 


K  (?)  -  1 


FOR  ?  —  oc  (y  —  o) 


(91) 


The  boundary  conditions  Eqs.  (90)  and  (91)  enable 
Eq.  (87)  to  be  integrated  once  the  viscosity  func¬ 
tion  G  is  specified.  Because  G(K)  is  a  fairly  com¬ 
plex  function  that  is  defined  by  an  implicit  relation 
(Eq.  86)  the  turbulent  energy  equation,  Eq.  (86), 
must  in  general  be  integrated  numerically.  A  simple 
exact  analytic  solution  satisfying  the  boundary 
conditions  can  be  found  for  the  special  case  that 
G  (K)  is  a  constant  equal  to  its  value  at  the  wall.  In 
this  case,  G(K)  -  1  and  K  (  )  is  found  by 
integration . 


K  (?)  3 


r  1  +  e-'^l 


(92) 


This  turns  out  to  a  good  approximation  even  when 
G  is  variable,  as  given  by  the  expression  in 
I.juboja  and  Rodila  and  serves  as  a  good  initial  es¬ 
timate  for  the  iterative  solution  for  K  (  -  )  in  the 
general  case. 


The  solution  for  K  ( ,-.)  determines  the  function 
B  (K).  appearing  in  the  differential  equation  for 


ii 


m 


E  (c).  This  equation  is  a  second-order  linear 
ordinary  differential  equation  which  can  be  solved  as 
initial  value  problem  using  far  field  boundary  con¬ 
ditions  determined  from  the  matching  conditions  to 
the  diffusion  layer  solution  for  c-»-0  (y-*-  <*> )  given 
by  Eq.  (79b).  Equation  (88a)  has  two  power  law 
type  solutions  for  *-0,  one  decaying  and  the 
other  unbounded.  In  order  to  match  to  the  dif¬ 
fusion  layer  the  unbounded  component  must  be  ex¬ 
cluded  by  imposing  the  far  boundary  condition , 


E(f )  —  0 


0  (y  —  »0 


The  solution  can  be  completely  determined  by  pre¬ 
scribing  the  coefficient  of  the  decaying  component. 
Because  of  the  arbitrary  coefficient,  e,  introduced 
in  the  definition  of  E  by  Eq  (84)  we  may,  without 
loss  in  generality ,  set  this  coefficient  to  unity  . 

Thus  the  second  boundary  condition  on  E(;)  is 
taken  as, 

lini  (E  fa*)  1 

e-o  (94) 

where  the  exponent,  a^,  is  determined  from  the 
differential  equation , 


where  a,  b  are  the  parameters  appearing  in  the 
outer  solution  as  defined  by  Eq.  (45).  These  con¬ 
ditions  enable  us  to  compute  E (  c )  as  a  simple  ini¬ 
tial  value  problem  by  marching  the  solution  from 
5  =  o  to  <*>  .  However,  we  found  it  more  conve¬ 
nient  to  solve  the  E  equation  as  a  two  point  bound¬ 
ary  value  problem  using  the  same  computer  code 
used  to  compute  K  (  O  .  The  velocity  profile 
u(y)  is  determined  implicitly  via  Eq.  (82)  and  the 
transformation  function ,  Y  ( <;) ,  is  evaluated  by 
a  numerical  integration  of  Eq.  (89).  The  scale 
function,  e,  for  the  dissipation  is  evaluated  by 
matching  the  transition  layer  solution  for  to  the 
far  field  limit  given  in  Eq.  (79b).  This  completes 
the  determination  of  the  solution  in  the  transition 
layer  in  terms  of  the  local  value  of  the  skin  friction . 
The  skin  friction  is  then  evaluated  by  carrying  out 
the  matching  of  transition  and  wall  layer  solutions. 

Limiting  Behavior  and  Evaluation  of  the  Skin 
F  riction 

The  parameter,  e,  scaling  the  dissipation  is 
determined  by  matching  the  transition  solution 
to  the  diffusion  layer.  For  C  —►0,  (y-*-« )  the  solu¬ 
tion  behaves  as, 

k  -  (6kw/G„)  r2  FOR  £  -  0  <y  —  <*)  (96) 


f  —  e£2 


FOR  £  —  0  <y  — «)  07) 


y  — (rwe"’  d'3*^2*1)  FOR£-0(y~°")  (98) 

where  G„  is  the  limiting  value  of  G  in  the  far  field 
and  is  given  by 

G *  -  GU()/GU„  (99a) 

where  Guw  Gu(P„/f„,  1)  (99b) 

and  d  is  the  constant,  =  Gu(0, 1)  (99c) 


[  18aG?/jl- 

d  ttCJ 


Equations  (96)  (98)  can  be  inverted  to  give  explic 

it  expressions  for  the  far  field  behavior  of  the  trun_ 
sition  solution  in  terms  of  the  physical  coordinate,  y. 
The  resulting  expressions  match  the  limiting  forms 
of  the  diffusion  layer  solution,  given  in  Eqs.  (79), 
if  the  arbitrary  function,  e,  scaling  the  dissipation 
is  given  the  value, 

e  =  df;b/a  B2/3l(x)  (]01) 

where  B  (x)  is  the  known  function  that  arises  in 
the  inner  limit  of  the  outer  solution  (Eq.  42). 

The  matching  to  the  diffusion  layer  solution 
and  the  boundary  condition  k-*-kw  completely  deter 
mines  the  solution  in  the  transition  layer  in  terms 
of  the  unknown  value  of  the  skin  friction,  tw.  This 
is  determined,  and  the  solution  completed  by  match¬ 
ing  the  transition  layer  solution  to  the  logarithmic 
law  of  the  wall. 

In  the  present  notation  the  law  of  the  wall  so¬ 
lution  in  the  inner  layer  is  expressed  in  the  form, 

U */u*  =  iq  F*  (yj ,  y*  =  yiq  Re,  (102) 

For  y+-»-  «  this  reduces  to  the  standard  log  law, 

U*/Uj*  -  iq(x)  [*•'  In  (yuT(x)  Re,)  +  CJ]  FOR  y*  -  *  ( 103) 

where  C+ is  the  standard  constant  (C+~5)  in  the 
logarithmic  law  of  the  wall.  For  „-►«>  (y-*-0)  the 
limiting  solution  for  the  turbulent  energy  approaches 
its  value  at  the  wall,  exponentially  in 


k„  -  k*[l  4  Q(e*^c)l 


FOR  £  -  «°  (104) 


while  the  dissipation  <.  and  coordinate  transforma 
tion ,  y  (  0  ,  approach  the  limits 


—  e  E.e-8*' 


-2*  /J*_  e-»-' 

e  E,  J  ck  0. 


FOR  £  -  ■»  (105) 


FOR  £  -  «> 


where  bw  is  a  constant  equal  to  the  limiting  value  of 
8  (K)  at  the  wall  and  is  given  by  the  following  ex¬ 
pression  involving  only  the  turbulence  model 
constants, 

Ar  =  (<vV<ce2  -  V  (107) 

The  parameter  Eg  is  a  constant  to  be  determined 
from  the  numerical  solution  for  E( :  )  for  — *•  ». 
Inverting  these  expressions  we  obtain, 

7~t3J2kJ'1  (108) 


u— *  — *-  [In (y e  Esk{ /tJ/2)  *  0(1)1 
K  i 

where  kj  must  satisfy  the  relation 


X2,  =  -  CM>  (110) 

This  is  just  one  of  the  standard  relations,  written  in 
the  present  notation,  used  to  fix  the  model  constants 
in  the  k  -  c  model  to  assure  that  it  satisfy  the  stan¬ 
dard  log  law  at  the  wall  (eg  see  Ref.  20).  Equations 
(  104),  (  108)  -  (  110)  demonstrates  that  the  transi¬ 
tion  layer  solution  matches  to  the  wall  layer  solution. 
Using  Eq.  (109)  the  transition  solution  for  the  ve- 
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locity  profile  near  the  wall  can  be  written  in  the  form , 

U  =  u,(x)  +7 ln(y dE,B2/3aK1  q(U1)/2*t*1  /a]  till) 

comparing  this  expression  to  the  outer  limit  of  the 
wall  layer  solution  given  by  Eq.  (103)  yields  the 
following  implicit  relation  for  the  skin  friction 

A'1  -  (1  +  2a'1)  In  A  =  (1  +  2a*1)lnuw(x)  -  In  9(x) 

^  [lna<1*))/4Re,  *  id  (112) 

where  the  skin  friction  is  given  in  terms  of  A  and 
uw  by 

rw/ui(x)  =  o)/2k2A:(x)  (113) 

and  9  (x)  is  the  function  defined  by, 

9<x)  =  dE,x;,2/aB2/3a(x)  (i14) 


In  the  far  field  limit,  x- 
following  limit 


'  ,  a  (x)  approaches  the 


0(x)  =  escj/ax*<,<1*a>/a  (115a) 


e„  =  dxi2/aEaB2/3a  (115b) 

and  Bg  is  the  constant  appearing  in  the  similarity 
expression  for  B(x)  given  in  Eq.  (60a).  Using  this 
relation  and  the  similarity  form  for  uw(x)  given  in 
Eq.  (31)  we  can  write  the  far  field  asymptotic  solu¬ 
tion  for  A(x)  in  the  form, 

A*1  -  (1  *  2a'1)  In  A  = 

Un(o(u,>/4RejCi/2e-ix(i-i)/2)  ,  kC*j  (116) 

where  C j  is  the  normalized  momentum  flux  defined  in 
Eq.  (32).  Equation  (116)  expresses  the  solution  for 
the  skin  friction  in  the  far  field  similarity  region  in 
terms  of  the  momentum  flux  constant ,  Cj ,  the  jet 
Reynolds  number.  Rej.  turbulence  model  constants 
and  two  constants  Bs  and  Es  which  are  determined 
by  numerical  solution  for  the  outer  (free  jet)  region 
and  transition  layer  equations.  It  is  interesting  to 
note  that  the  local  skin  friction  coefficient,  cp  =  tw/ 
(1/2)  uw9,  is  exactly  constant  for  a  radial  wall  jet 
(j  =  1)  and  varies  logarithmically  in  x  for  the  planar 
jet  (j  =  0) .  The  above  formula  for  the  skin  friction 
admits  an  interesting  similarity  that  will  be  dis¬ 
cussed  in  the  next  section. 


The  similarity  form  for  the  radial  and  planar 
outer  layer  and  the  transition  layer  have  been  solved 
numerically  for  Gy  =  1  as  well  as  for  the  algebraic- 
stress  model  of  Ljuboja  and  Rodi^,  here  denoted  by 
G  #  1.  Product  composites  of  the  form 
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where  CP  denotes  the  common  part  of  the  subscript¬ 
ed  layers,  are  used  to  compare  the  present  asymp¬ 
totic  solutions  with  respect  to  the  numerical  solutions 
of  bjuboja  and  Rodi.13  These  results  are  also  com¬ 
pared  with  "universal"  correlations  of  experimental 
data  proposed  by  Narasimba  et  al.  for  plane  wall 
jets. 

The  equations  governing  the  outer  layer  are 


solved  by  the  method  of  Paullay  et  al.19  Application 
of  an  independent  variable  transformation,  5  =  5(n) 
such  that 

dn  =  G  — 

d£  “  11  (118) 

acts  to  decouple  the  kinetic  energy  and  dissipation 
equations  ns  well  as  to  reduce  much  of  the  non¬ 
linearity  of  the  system.  In  the  new  coordinate  the 
bounds  +  »  replace  a  =0,  riband,  the 

location  of  the  finite  edge  is  found  a  posteriori 
by  quadrature,  viz. 


r  “  r-2 

’).  =  g“TT  d* 


The  numerical  solution  is  obtained  by  central  differ¬ 
encing  the  transformed  equations  and  using  the  edge 
and  wall  asymptotic  behavior  as  two-point  boundary 
values.  Repeated  application  of  tridiagonal  solvers 
to  the  system  of  equations  allows  for  their  conver¬ 
gence  to  levels  of  local  residuals  no  greater  than 

10-11. 

Results  for  the  outer  solution  in  the  planar 
case  with  Gu  =  1  (Fig.  3)  are  compared  with  the 
corresponding  free  jet  solutions.!9  For  reasons 
discussed  previously ,  the  radial  and  planar  wall  jet 
solutions  have  the  identical  similarity  form  as  the  free 
jet  flow.  It  is  apparent  that  the  effect  of  the  wall 
is  to  reduce  the  velocity  half-width  by  about  30®. 

The  behavior  of  the  mean  and  turbulent  variables 
are  essentially  the  same  in  the  edge  region  of  the 
layer  but  differ  significantly  near  rv^O  where  wall 
jet  dissipation  becomes  unbounded  and  the  turbulent 
energy  and  eddy  viscosity  go  to  zero.  The  results 
for  radial  jets  are  qualitatively  the  same  as  for  the 
planar  case. 

The  mean  velocity  deficit,  u/um,  profile  in  the 
diffusion  layer  (Fig.  4)  is  given  analytically  by 
Eq  ( 57)  and  it  exhibits  the  local  maximum  that 
characterizes  wall  jets. 

Solutions  to  the  transition  layer  must,  in  gen¬ 
eral,  be  obtained  numerically  and  again  central  dif¬ 
ferences  are  employed  to  represent  Eqs  (87,88) 
subject  to  boundary  conditions  Eq  (90,91)  and  initial 
conditions  Eq  (93,94),  respectively.  In  practice  K"1 
is  used  as  the  dependent  variable  in  Eq  (87)  and  the 
difference  equations  are  solved  from  5=0  (y-+»)  to 
some  finite  cutoff,  5  =  5c>  chosen  large  enough  for 
a  two  term  large  5  expansion  to  prevail.  Since  the 
dissipation  equation,  Eq  (88),  is  homogeneous,  it 
may  be  solved  as  a  two-point  boundary  value  problem 
with  E  set  to  zero  at  5  =  0.  Eq  (93)  and  to  an  arbi¬ 
trarily  prescribed  value  at  the  edge  c  =  5c .  The  re¬ 
sulting  numerical  solution  is  then  rescaled  to  comply 
with  Eq  (94).  This  procedure  enabled  us  to  solve 
for  K  and  E  simultaneously  on  the  same  mesh  to  the 
same  level  of  precision.  The  numerical  results  are 
indistinguishable  from  the  analytical  solution  for 
K  (5).  Eq  (92).  that  is  available  for  Gu  =  1.  The 
more  general  algebraic  stress  model!9  (i.e.  Gu  #  1) 
gives  results  (Fig.  5)  that  are  very  close  to  those 
obtained  with  the  standard  (i.e.  G,j  =  1)  turbulence 
model.  It  should  be  recalled  that  the  Y(  .)  curve  is 
actually  a  stretched  velocity  profile  (Eqs  82,89). 

The  numerical  solutions  for  the  outer  layer  and 
the  transition  layer  provide  the  constants  (Table  1) 
required  for  asymptotic  matching  with  the  diffusion 
layer.  It  has  been  established  that  all  numerical 
solutions  exhibit  the  correct  asymptotic  behaviors  for 
■  -  rc,  0  and  5  *  0,  «■  with  the  exponents  pro 

dieted  to  within  1®  of  their  analytical  values. 
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TABLE  1  Computed  Constants  in  the  Similarity  Solutions 


TURBULENCE 

MODEL 

GEOMETRY 

B. 

E, 

'2 

STANDARD  WITH  PLANAR 

.0635 

6.30x1 0*3 

18.96 

.1899 

-  0.09 

RADIAL 

.1106 

1.218x10'2 

18.96 

.1308 

MOOIFIED  WITH 

PLANAR 

.0337 

1.370x10"3 

12.45 

.0932 

CM  -CM  (p/e .  f) 

! 

RADIAL 

.0575 

3.05x10‘3 

12.45 

.0591 

In  order  to  usi  these  results  to  generate  the 
solution  profiles  for  specific  wall  set  flows  the  skin 
friction  must  first  be  evaluated  from  Eq  (116).  It  is 
convenient  to  rewrite  this  equation  in  the  form 

A*1  -  (1  +2a‘l)  In  A  =  i  In  ,3  +  {iln(al2e;2)  +  xC*}  (120) 

WherC  ^  =  =  (^j»  <«*1/2x»U1  ( 12D 

This  formulation  demonstrates  that  for  the  radial 
wall  jet  the  skin  friction  coefficient ,  A  ,  is  a  con¬ 
stant  that  depends  only  on  the  jet  Reynolds  number 
and  turbulence  model  constants;  for  a  plane  wall  jet  A 
is  variable  that  depends  on  the  logarithm  of  the 
streamwise  coordinate,  x.  The  solution  to  Eq  (120) 
is  presented  (Fig.  6)  for  radial  and  plane  wall  jets 
and  for  GU  =  1,GU*1.  As  expected,  A  decays 
slowly  with  J.  For  the  data  to  be  examined  here  J 
is  in  the  109  <  J  <  10 12  range  where  A  =  0.1. 

13 

A  solution  corresponding  to  Ljuboja  and  Rodi's 
numerical  result  at  r*/b*  =  200  (x  =  60)  (Rej  =  1.8  x 
104,  J  =  6.48  x  1010)  is  developed  using  their  alge¬ 
braic  Reynolds  stress  model  (Gu*  1).  The  composite 
solution  for  velocity  is  presented  in  Fig.  7,  along 
with  the  outer  and  transition  layer  calculations  in 
order  to  indicate  the  overall  nature  of  the  matching. 
The  large  gap  between  the  outer  and  composite 
solutions  apparent  in  this  result  is  an  indication  of 
the  importance  of  higher  order  terms  in  the  outer 
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Fig.  7  Structure  of  Composite  Solution  -  Plane  Wall  Jet 
(G^  *  1,  A  =  0.100) 


solution.  The  velocity  profile  merges  with  the  law  of 
the  wall  behavior  at  about  y  +  -  20  (Fig.  8)  which  is 
close  to  the  patching  location  used  by  Ljuboja  and 
Rodi.  Since  the  law  of  the  wall  region  extends  out  to 
about  y  +  *20  these  results  indicate  that  typical  solu¬ 
tions  exhibit  a  logarithmic  law  of  the  wall  behavior 
only  over  a  very  limited  interval.  This  marginal  log 
law  region  is  consistent  with  the  experimental  ob¬ 
servations8.  10 , 12  discussed  previously.  The  mean 
velocity,  turbulent  shear  stress  and  energy  profiles 
obtained  13  by  a  parabolic  marching  integration  of  the 
partial  differential  equations  are  in  reasonable  agree¬ 
ment  with  composite  for  field  solution  calculated  here 
(Fig.  9). 
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F  ig.  6  Friction  Velocity  a«  a  Function  of  Wall  Jet  Momentum 
Parameter  from  Transition  Layer  Matching  with 
Law  of  the  Wall. 


Fig.  8  Composite  Solution  Match  with  Law  of  the  Wall  -  Plane 
Wall  Jet  (G^  =  1,  A  «  0.100,  k  *  0.43,  C+  =  5.0) 
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Fig.  9  Comparison  of  Asymptotic  Composite  Solution  with 
Parabolic  Matching  Result  for  Plane  Wall  Jet 

Agreement  is  good  in  the  region  near  the  wall 
and  deteriorates  toward  the  edge  of  the  wall  jet,  which 
is  another  indication  of  the  need  for  including  higher 
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order  terms  in  the  outer  layer.  The  model  constants 
used  in  the  algebraic  stress  model  of  Ref.  13  are  not 
consistent  with  the  law  of  wall  compatability  condition 
given  by  Eq  (110).  In  order  to  have  the  outer  layer 
turbulence  model  conform  to  that  of  Ljuboja  and 
Rodi,  Eq  (110)  is  used  throughout  to  locally  deter¬ 
mine  ce2  as  a  function  of  G^.  Numerical  experi¬ 
ments  carried  out  here  and  by  Rodi21  indicate  neg¬ 
ligible  effects  due  to  this  slight  inconsistency  in 
their  model. 

A  dimensional  analysis  performed  by  Narasimha 
et  al.  18  suggests  that  plane  wall  jets  are  described 
by  their  initial  momentum  and  the  kinematic  viscosity . 
Together  these  quantities  provide  unique  length  and 
velocity  scales  [L*]  =  [y*2/j*],  [u*]  =  [J*/y*1- 
For  the  radial  wall  jet  this  is  not  the  case  since  J */ 
Y*2  is  itself  dimensionless.  The  results  obtained 
here  conform  to  the  analysis  of  Narasimha  et  al.  The 
plane  wall  jet  shear  stress  is  determined  only  by  the 
value  of  ,0  (Eq  120,  Fig  6).  The  streamwise  behavior 
of  the  scaled  mean  and  turbulence  quantities  are 
universal  functions  of  ))  which  serves  as  a  scaled 
streamwise  coordinate  in  the  planar  case.  (Fig. 
lOa-d).  Comparison  of  the  present  planar  results 
with  the  proposed  universal  curves*®  show  good 
agreement  for  velocity  maximum  decay  (Fig.  10a). 

The  present  results  for  the  maximum  velocity  were 
computed  from  the  composite  solution  given  in  Eq 
(117)  rather  than  from  the  analytic  expression  given 
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by  the  diffusion  layer  solution,  Eq  (64).  The  velo¬ 
city  half-widths  are  underpredicted  by  30%  (Fig. 

10b).  The  calculation  of  Ljuboja  and  Rodi  agree 
with  Narasimha  giving  further  cause  to  surmise  that 
higher  order  term  terms  in  the  outer  layer  are 
.mportant.  The  velocity  half-width  is  expected  to  be 
quite  sensitive  to  such  terms.  The  height  above  the 
wall,  y*  ,  of  the  velocity  maximum  is  considerably 
underpredicted  by  the  composite  solution  (Fig.  10c) 
and  this  is  partly  explained  by  the  use  of  an  eddy 
viscosity  model  which  forces  this  point  to  be  collo¬ 
cated  with  the  zero  shear  stress  point.  Since  the 
vanishing  shear  stress  location  is  observed  to  fall 
between  the  wall  and  the  velocity  maximum ,  some 
underprediction  is  expected.  Still,  the  computations 
of  Ljuboja  and  Rodi  are  considerably  closer  to  the 
universal  curve  of  Narasimha  et  al.  The  behavior  of 
the  wall  shear  stress  (Fig.  lOd)  is  of  particular  note. 
Narasimha  et  al  give  no  such  result  but  Hammond9 
correlates  wall  shear  stress  with  a  scaling  based  upon 
slot  height.  Here  the  Hammond  result  for  skin  friction 
is  rescaled  and  compared  with  our  analytic  solution , 
Eq(120),  the  calculations  of  Ljuboja  and  Rodi1-5  and 
the  observations  of  Bradshaw  and  Gee22.  The 
agreement  of  the  latter  with  the  proposed  curve  of 
Hammond  is  a  direct  result  of  our  rescaling.  The 
Bradshaw  and  Gee*-  result  did  not  correlate  well  with 
the  original  scaling  due  to  Hammond9.  The  present 
solution  shows  the  correct  variation  of  the  wall  shear 
with  distance  but  overpredicts  its  level  by  about  40% 
again  suggesting  a  sensitivity  of  Eq  (120)  to  higher 
order  terms  in  the  outer  layer.  Nevertheless,  the 
results  presented  here  confirm  (at  least  to  lowest 
order)  that  the  scaling  of  Narasimha  et  al. 18  does 
indeed  provide  a  universal  correlation  for  the  far  field 
of  a  plane  wall  jet.  Interestingly,  radial  wall  jets  have 
already  been  shown  to  possess  far  field  composite 
profiles  (i.e.  across  the  entire  wall  jet)  that  depend 
only  upon  their  initial  momentum  (i.e.  idependent  of 
streamwise  distance).  Although  plane  jet  far  field 
profiles  vary  slowly  with  streamwise  distance,  their 
streamwise  decay  is  described  by  a  universal  function 
of  a  single  variable  dependent  on  initial  momentum 
and  distance  while  those  of  radial  wall  jets  are  not. 
The  lack  of  such  universality  for  radial  wall  jets  as 
well  as  the  spareseness  of  far  field  observations  pre¬ 
cludes  their  comparison  with  the  composite  solution 
at  the  present  time. 

VI .  Conclusions  and  Discussion 

The  present  paper  describes  a  new  analytical 
approach  for  analyzing  turbulent  plane  and  radial 
wall  jets.  The  method  is  based  on  a  systematic  two 
parameter  asymptotic  expansion  using  a  k-t  turbu¬ 
lence  model.  One  of  the  parameters,  y,  is  related 
to  the  friction  velocity  and  controls  the  Reynolds 
number  effects  induced  by  the  solid  wall;  the  other 
parameter,  a,  is  equal  to  the  model  constant  appear¬ 
ing  in  the  eddy  viscosity  formula  and  controls  the 
turbulent  shear  stress  levels  in  the  outer  free  jet 
part  of  the  wall  jet.  The  a  expansion  leads  to  the 
usual  boundary  layer  approximation  while  the  y 
expansion  leads  to  a  four  layer  structure  of  the 
wall  jet  for  y->0.  The  theory  is  valid  in  a  region 
downstream  from  initial  zone  where  the  wall  jet  is 
formed.  He  have  determined  the  equations  and 
matching  conditions  governing  the  leading  term 
of  the  expansion  in  each  layer.  The  theory  leads  to 
explicit  expressions  for  the  magnitude  and  position  of 
the  maximum  velocity  in  the  wall  jet  and  for  the  shear 
stress  at  the  wall.  Similarity  solutions  have  been 
developed  for  each  layer  to  describe  the  limiting 
behavior  of  the  wall  jet  far  downstream.  Complete 
solutions  of  the  leading  order  similarity  equations 


in  each  layer  have  been  obtained  for  plane  radial 
jets  with  both  the  standard  k-<_  turbulence  model 
(constant  cu)  and  for  the  Ljuboja  and  Rodi 
algebraic  reynolds  stress  model.  The  present 
solutions  for  the  plane  wall  jet  where  compared  with 
Ljuboja  and  Rodi's  numerical  solution  of  the  bound 
ary  layer  and  k-  c.  model  equations  and  with 
Narasimha,  Yenga  Narayan  and  Parthasarathy's 1 
experimental  correlations . 

The  main  conclusions  of  the  present  work  are; 

•  The  wall  jet  develops  a  four  layer  structure 
for  Rej-n”  ■  The  multi  layer  structure  is  a 
consequence  of  the  turbulent  shear  stress 
in  the  outer  layer  being  large  compared  to 
the  shear  stress  levels  in  the  inner  w;ill 
layer.  We  believe  this  structure  is  a  gen 
eral  feature  of  other  turbulent  shear  Hows 
that  have  relatively  high  turbulent  shear 
stress  levels  in  the  outer  region. 

•  The  outer  region  of  a  wall  jet  is  similar  to 
a  free  jetT  The  present  theory  shows  that 
the  overall  momentum  flux  in  a  wall  jet  is 
constant  to  lowest  order,  us  it  is  in  a  free 
jet  and  that  the  similarity  form  of  the  outer 
region  of  the  wall  jet  is  similar  to  a  free 

jet  with  the  same  velocity  and  length  scales 
to  lowest  order.  The  present  theory  pro 
vides  a  theoreticiil  explanation  of  why  the 
observed  variation  of  the  velocity  and 
length  scales  of  wall  jets  are  so  close  to 
those  of  free  jets.  Nevertheless,  the  pres¬ 
ence  of  the  solid  surface  is  felt  in  the 
outer  flow  through  the  boundary  conditions 
on  the  turbulent  energy  and  dissipation . 
Thus  the  outer  solution  in  the  present  the¬ 
ory  is  not  a  free  jet  as  in  Glauert's  theory. 
The  behavior  or  the  outer  solution  near  a 
wall  is  dominated  by  the  vanishing  of  the 
turbulent  viscosity  at  the  wall  and  conse¬ 
quently  the  profiles  in  the  outer  solution 
differ  significantly  from  the  symmetrical 
profiles  of  a  free  jet . 

•  The  structural  form  of  the  wall  jet  solution 
is  sensitive  to  the  turbulence  model.  It  is 
the  use  of  the  dissipation  equation  that 
leads  to  the  requirement  for  a  four  layer 
structure.  The  length  scale  in  the  "extra" 
transition  region  is  sensitive  to  the  model 
constants  in  the  k-  e  model  through  their 
influence  on  the  parameter  a.  The  four  layer 
structure  for  this  type  of  flow  (i.e.  large 
turbulent  shear  stress  levels  in  outer  part 
of  shear  layer)  may  be  a  general  property 
of  the  dissipation  equation.  Although  not 
demonstrated  here,  it  can  be  shown  that 
for  a  prescribed  algebraic  length  scale  the 
extra  transition  region  is  not  required  and 
a  three  layer  structure  evolves.  The  con¬ 
sequences  of  this  observation  for  turbu¬ 
lence  modeling  remains  to  be  established. 

•  The  present  theory  indicates  dominance  of 
turbulent  diffusion  near  the  walh  The 
diffusion  layer  is  a  productionless  zone 
that  is  dominated  by  a  balance  of  diffusion 

.  and  dissipation  to  lowest  order.  Again, 

this  may  be  a  general  feature  of  flows  with 
relatively  large  turbulent  shear  stress  lev¬ 
els  in  the  outer  region  . 


•  The  velocity  maximum  occurs  in  the  diffu¬ 
sion  layer.  Solution  profiles  in  the  diffu 
sion  layer  are  described  by  simple  analytic 
expressions  that  follow  from  the  linear 
shear  stress  variation  across  the  layer  and 
from  the  balance  of  diffusion  and  dissipa¬ 
tion.  Simple  expressions  are  found  for  the 
magnitude  and  position  of  the  maximum 
velocity  in  the  wall  jet.  The  expression 
for  the  maximum  velocity  is  exactly  similar 
for  a  radial  wall  jet  but  is  not  for  a  plane 
wall  jet  where  there  is  an  additional  weak 
logarithimic  dependence  on  the  streumwise 
distance. 

•  T  he  wall  jet  has  an  asymptotic  similarity 
form  that  describes  the  flow  far  down¬ 
stream.  The  leading  order  terms  are  ex¬ 
actly  similar  for  x->-“.  Deviations  from 
exact  similarity  arise  from  higher  order 
terms  in  the  asymptotic  expansions.  The 
present  theory  accounts  only  for  leading 
order  similar  terms. 

•  The  present  theory  is  consistent  with  a 
conventional  law  of  the  wall  in  the  inner¬ 
most  layer.  Because  of  the  rapid  changes 
across  the  diffusion  and  transition  layers 
the  extent  of  the  overlap  region,  in  which 
the  velocity  profile  is  well  represented  by 
a  conventional  logarithmic  law  of  the  wall, 
is  small,  which  is  generally  consistent  with 
experimental  observations.  The  present 
theory  shows  that  this  diminished  overlap 
region  is  due  to  the  large  gradients  of 
shear  stress  normal  to  the  wall.  The  pres¬ 
ent  theory  indicates  that  this  reduces  the 
extent  of  the  log  law  region  but  does  not 
eliminate  it. 

•  The  theory  leads  to  a  solution  for  the  skin 
friction  in  the  form  of  a  standard  log  law. 

A  general  explicit  expression  for  the  skin 
friction  is  deduced  for  the  intermediate 
streamwise  zone  (i.e.  x  =  0(1))  of  the  wall 
jet  where  the  outer  solution  is  described 
by  partial  differential  equations.  The  so¬ 
lution  involves  one  unknown  function  B(x) 
that  must  be  determined  from  the  numerical 
solution  of  the  outer  layer  partial  differen¬ 
tial  equations.  A  conipletely  determined 
solution  for  the  skin  friction  is  deduced  for 
the  asymptotic  far  field  similarity  region 

(x  -*•  «).  The  loc(d  skin  friction  coefficient, 
defined  in  terms  of  local  values  of  Uw  (x) 
or  Um  (x) ,  is  shown  to  be  constant  for 
radial  jets  and  to  depend  on  log  x  for  plane 
jets. 

•  The  present  theory  provides  the  basis  for 
the  scaling  analysis  that  Narasimha,  Yenga 
Narayan,  and  Parthasarathy~t used  to 
successfully  correlate  plane  wall  jet  data. 
The  present  theory  leads  to  the  same  sim¬ 
ilarity  variables  deduced  by  these  authors 
on  the  basis  of  dimensional  analysis  and 
shows  that  the  plane  wall  jet  solution  de 
pends  only  on  a  suitably  defined  non- 
dimensional  momentum  flux.  The  theory 
explains  why  a  similar  type  dimensional 
analysis  does  not  work  for  radial  wall 
jets. 


•  Comparisons  of  the  present  solutions  for 
plane  w;dl  jet  profiles  shows  relatively 
good  agreement  with  l.juboja  and  Rodi's 
numerical  solutions  of  the  boundary  layer 
equations.  The  solution  for  the  decay  of 
the  maximum  velocity  also  shows  good 
agreement  with  Ljuboja  and  Rodi's  compu 
tations  and  with  the  experimental  correla¬ 
tions  of  Ref.  18. 

The  present  paper  was  limited  to  the  develop 
ment  of  the  overall  asymptotic  approach  and  to  the 
computation  of  the  leading  order  terms.  Initial 
comparisons  of  the  present  theory  with  "exact" 
numerical  calculations  ^  and  with  experimental 
correlations  ^  for  plane  wall  jets  shows  that  it 
correctly  predicts  many  features  of  these  flows. 
However,  these  comparisons  also  indicate,  (I’igs. 

10b,  c) ,  that  the  present  first  order  theory  signif 
icantly  underpredicts  the  spread  rates  61/2  and 
0  u]  of  both  the  half  velocity  and  maximum  velocity 
and  maximum  velocity  points,  respectively.  A1 
though  the  overall  variation  with  downstream  dis¬ 
tance  is  correctly  predicted,  due  mainly  to  sim¬ 
ilarity  considerations,  the  absolute  magnitudes  are 
not.  The  same  is  true,  to  a  lesser  degree,  for  the 
skin  friction  predictions  (Fig  lOd).  These  discrep¬ 
ancies  are  indicative  of  the  importance  of  the  ne¬ 
glected  second  order  terms  in  the  expansions.  The 
relatively  large  gap  between  the  outer  and  composite 
solutions  for  the  velocity  apparent  in  the  results  of 
Fig.  7  is  indicative  of  the  magnitude  of  neglected 
second -order  effects. 

The  present  theory  is  based  on  the  k-e  turbu 
lence  model,  which  employs  an  eddy  viscosity  model 
for  the  turbulent  shear  stress.  Although  the  k-e 
model  cannot  adequately  describe  the  details  of  the 
flow  near  the  points  of  maximum  velocity  and  zero 
shear  stress,  Ljuboja  and  Rodi  show  that  their 
k-t  formulation  can  predict  most  other  features  of 
wall  jet  flows  on  flat  surfaces,  where  these  two 
points  are  close  together.  As  indicated  by  Gibson 
and  Younisl4,  the  points  of  maximum  velocity  and 
zero  shear  stress  become  more  widely  separated  on 
curved  surfaces;  these  flows  cannot  be  well  pre¬ 
dicted  with  an  eddy  viscosity  formulation.  How¬ 
ever,  Gibson  and  Younis  show  that  a  significantly 
improved  description  of  the  How  near  a  velocity 
maximum  can  be  achieved  with  a  Reynolds  stress 
model  (RSM)  of  turbulence  that  incorporates  a  wall 
damping  term  of  the  type  used  by  Ljuboja  and  Rodi. 
The  asymptotic  approach  described  in  the  present 
paper  can  be  directly  applied  to  the  Gibson  and 
Younis  RSM  equations.  Since  the  RSM  employs  the 
same  dissipation  equation  used  in  the  present  study 
the  theory  will  lead  to  the  same  four  layer  structure 
as  in  present  work.  The  use  of  a  RSM  can  be  ex¬ 
pected  to  give  a  significantly  improved  description 
of  the  solution  in  the  diffusion  layer.  It  would  be 
useful  to  determine  whether  the  four  layer  struc¬ 
ture  is  a  general  feature  of  flows  with  an  outer 
region  of  large  Reynolds  shear  stress  or  if  it  is 
simply  a  consequence  of  the  particular  form  chosen 
for  the  dissipation  equation. 

We  have  carried  out  only  limited  comparisons 
of  the  theory  with  experimental  data.  It  is  certain¬ 
ly  desirable  to  continue  these  studies  and  to  per 
form  more  extensive  comparisons.  Howver,  before 
proceeding  with  futher  applications,  it  would  appear 
best  to  compute  the  second -order  terms  to  determine 
if  they  improve  the  predictions  of  spread  rate  and 
skin  friction  levels. 


18 


& 


«  V*  ~Tr 


Ufl 

**> 


Acknowledgement  11. 

The  research  reported  in  this  paper  was 
partially  supported  by  the  United  States  Office  of 
Naval  Research  under  contract  N00014-81-C-0549  12. 


VIII.  References 

1.  Townsend  A. A.  ■  The  Structure  of  Turbulent 
Shear  Flow .  2nd  ed.  ,  Cambridge  University 
Press,  New  York,  1976,  pp.  268-271. 


!:.i  n"  u 1 


»  !7  .  r;  ■.-'  V  W.' 


Meyers,  G.E.,  Schauer,  J.J.  and  Eustis,  R.W., 
"The  Plane  Turbulent  Wall  Jet,  part  I,"  Trans. 
ASME,  J.  Basic  Eng.,  85,  1963,  pp.  47-54. 

Bradshaw,  P.  and  Gee,  M.T.,  "Turbulent  Wall 
Jets  with  and  without  an  External  Stream," 
ARCR&M  3252,  1960. 

Ljuboja,  M.  and  Rodi,  W.,  "Calculation  of 
Turbulent  Wall  Jets  with  an  Algebraic  Reynolds 
Stress  Model,"  Trans.  ASME,  J.  Fluid  Eng., 

102.  1980,  pp.  350-356. 


US 


A  «L 


0 


■- .' 


2.  Ibid,  pp.  195-230,  pp.  262-276. 

3.  Rotta,  J.C.,  "Turbulent  Boundary  Layers  in 
Incompressible  Flow,"  Progress  in  Aero.  Sci, 
2ed.,  A.  Ferri,  D.  Kuchemann  and  L.H.G. 
Sterne,  Pergamon  Press,  Oxford,  1962,  pp 
1-219. 

4.  Rodi,  W.,  The  Prediction  of  Free  Turbulent 
Boundary  Layers  by  Use  of  a  Two-Equation 
Model  of  Turbulence,  Ph.D.  Thesis,  Univer¬ 
sity  of  London,  1972. 

5.  Schlichting,  H.,  Boundary  Layer  Theory,  7th 
ed.,  McGraw  Hill  Book  Co.,  New  York  1979, 
Chapter  XXIV. 

6.  Glauert,  M.B.  "The  Wall-Jet"  J .  Fluid  Mech .  , 

_1_,  1956,  pp.  625-643. 

7.  Bakke,  P.  ,  "An  Experimental  Investigation  of 
a  Wall  Jet,"  J.  Fluid  l.Iech.  ,  2,  1957,  pp. 467- 
472. 

8.  Poreh,  M.,  Tsuei,  Y.G.  and  Cermak,  J.E., 

"Investigation  of  a  Turbulent  Radial  Wall  Jet," 
Trans.  ASME,  J.  Appl.  Mech.,  89,  1967,  pp. 
437-463. -  - 

9.  Hammond,  G.P.,  "Complete  Velocity  Profile 
and  "Optimum"  Skin  Friction  Formulas  for  the 
Plane  Wall-Jet,"  Trans.  ASME,  J.  Fluids  Eng., 

104,  1982,  pp.  5jnnr 

10.  Schwarz,  W.H.  and  Cosart,  W.P.,  "The  Two 
Dimensional  Turbulent  Wall  Jet,"  J.  Fluid  Mech, 
U>,  1961,  pp.  481-495. 


14.  Gibson,  M.M.  and  Younis,  B.A.,  "Modeling  the 
Curved  Turbulent  Wall  Jet,"  A1AA  J ■ ,  20, 

1982,  pp.  1707-  1712. 

15.  Launder,  B.E.  and  Rodi,  W.,  "Turbulent  Wall 
Jet,"  1980-81  AFOSR-HTTM- Stanford  Confer¬ 
ence  on  Complex  Turbulent  Flows,  ed.  S.J. 
Kline,  B.J.  Cantwell  and  G.M.  Lilley ,  Stanford 
University,  1981,  pp.  434-456. 

16.  Adamson,  T.C.  and  Messiter,  A.F.,  "Analysis 
of  Two-Dimensional  Interactions  Between  Shock 
Waves  and  Boundary  Layer,"  Ann  Rev.  Fluid 
Mech ,  12,  1980,  pp.  103-138. 

17.  Melnik,  R.E.,  "Turbulent  Interactions  on  Air¬ 
foils  at  Transonic  Speeds  -  Recent  Develop¬ 
ments,"  AGARD  Conference  of  Computation  of 
Viscous  Inviscid  Interactions,  CP-291,  Paper 
No.  10,  Colorado  Springs,  Co,  Sept.  1980. 

18.  Narasimha,  R.,  Yenga  Narayan,  K.  and 
Parthasarathy ,  S.P.  ,  "Parametric  Analysis  of 
Turbulent  Wall  Jets  in  Still  Air.  ."Aeronautical 
J_. ,  77,  1973,  pp.  355-359. 

19.  Paullay,  A.J.,  Melnik,  R.E.,  Rubel,  A., 
Rudman,  S  and  Siclari,  M.J.,  "Similarity 
Solutions  for  Plane  and  Radial  Jets  Using  a 
K-E  Turbulence  Model"  (submitted  to  JFM). 

20.  Rodi,  W.  ,  Turbulence  Models  and  Their  Appli¬ 
cation  in  Hydraulics,  IAHR-Section  on  Funda¬ 
mentals  of  Divison  II ;  Experimental  and  Math¬ 
ematical  Fluid  Dynamics,  1980,  p  28. 

21.  Rodi,  W.,  Private  Communication  Nov.,  1982. 


0 1  ST IRBUT ION  LIST  FOR  UNC L ASS  I c I E n 
TECHNICAL  REPORTS  AND  REPRINTS  ISSUED  UNDER 
CONTRACT  N00014-81-C-0549  TASK  061-251 

All  addresses  receive  one  copy  unless  otherwise  specified. 


Technical  Library 
8ui Idioq  313 

Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  MD  21005 

Mr.  Aviars  Cel  mi  ns 
Ballistic  Research  Laboratory 
Ballistic  Modelling  Di v i d i on 
Aberdeen  Proving  Ground,  MO  21005 

Dr.  P.J.  Roache 

Ecodynami cs  Research  Associates, 

Inc. 

D.0.  Box  8172 
Albuquerque,  I'M  87108 

Defense  Technical  Information  Center 
CAmeron  Station,  Building  5 
Alexandria,  VA  22.114  12  copies 

Library 
Naval  Academy 
Annapolis,  MD  21402 

Di rector.  Tactical  Technology  Office 
Defense  Advanced  Research  Projects 
Agency 

1400  Wi 1  son  Boul evard 
Arlington,  VA  22209 

Code  200B 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217 

Code  438 

Office  of  Naval  Research 

800  N.  Quincy  Street 

Arlington,  VA  22217  2  copies 

Dr.  J.L.  Potter 

Deputy  Director,  Technology 

von  Karman  Gas  Dynamics  Facility 

Arnold  Air  Force  Station,  TN  37389 


Professor  J.  C.  Wu 
School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  GA  3033? 

Li  brary 

Aerojet-General  Corporation 
6352  North  Irwindale  Avenue 
Azusa,  CA  91702 

NASA  Scientific  and  Technical 
Information  Faci 1 ity 
P.0.  Box  8757 

Bal timore/Washiogton  International 
Airport,  MD  21240 

Dr.  K.  C.  Wang 
College  of  Engineering 
San  Diego  State  University 
San  Diego,  CA  92182 

Professor  A..J.  Chorin 
Department  Mathematics 
University  of  California 
Berkeley,  CA  94720 

Professor  M.  Hoi t 

Department  of  Mechanical  Engineering 
University  of  California 
Berkeley,  CA  94720 

Dr.  H.  R.  Chaplin 
Code  1600 

Oavid  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  MD  20084 

Dr.  Hans  Lugt 
Code  1802 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  MD  20084 


Or.  Francois  Frenkei 1 
Code  1802 

David  VI.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  Ml)  20084 

Dr.  T.C.  Tai 
Code  1806 

David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  MD  20084 

Professor  C.H.  Lewis 
Department  of  Aerospace  and  Ocean 
Fn  qi  nee1"  i  nq 

Virginia  Polytechnic  Institute  and 
State  University 
81 acksburg,  VA  24061 

Professor  A.H.  Nayfeh 
Department  of  Engineering  Science 
Virginia  Polytechnic  Institute  and 
State  University 
Blacksburg,  VA  24051 

Commanding  Officer 
Office  of  Naval  Research 
Eastern/Central  Regional  Office 
666  Summer  Street,  Bldg.  114, 
Section  0 
Boston,  MA  02210 

Dr.  J.C.  Erickson,  Jr. 

CALSPAN  Corporation 
Advanced  Technology  Center 
P.0.  Box  400 

Buffalo,  New  York  14225 

Dr.  T.J.  Falk 
CALSPAN  Corporation 
Advanced  Technology  Center 
P.0.  Box  400 

Buffalo,  New  York  14225 

Dr.  C.  Wit!  iff 
CALSPAN  Corporation 
Advanced  Technology  Cnter 
P.0.  Box  400 

Buffalo,  New  York  14225 


Professor  R.  F.  Probstein 
Department  of  Mechanical  Engineering 
Massachusetts  Institute  of 
Technol ogy 

Cambridge,  MA  02139 

Commanding  Officer 

Officer  of  Naval  Research  Branch 

Office 

536  South  Clark  Street 
Chicago,  IL  60605 

Code  753 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Mr.  J.  Marshall 
Code  4063 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Professor  R.T.  Davis 
Department  of  Aerospace  Engineering 
University  of  Cincinnati 
Cincinnati,  OH  45221 

Professor  S.G.  Ruoin 
Department  of  Aerospace  Engineering 
and  Appl ied  Mechanics 
University  of  Cincinnati 
Cincinnati,  OH  45221 

Library  MS  60-3 
NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 

Dr.  J.D.  Anderson,  Jr. 

Chairman,  Department  of  Aerospace 
En  gi neeri ng 

College  of  Engineering 
University  of  Maryland 
College  Park,  MD  10742 

Professor  0.  Burggraf 
Department  of  Aeronautical  and 
Astronautical  Engineering 
Ohio  State  University 
1314  Kinnear  Road 
Columbus,  OH  43212 


Technical  Library 
Naval  Surface  Wapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22443 

:)r.  F.  Moore 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22443 

Technical  Library  2-51131 
LTV  Aerospace  Corporation 
P.0.  Box  5907 
Dallas,  IX  75222 

Library,  United  Aircraft  Corporation 
Research  Laboratories 
Si  1 ver  Lane 

Fast  Hartford,  CT  06108 
Dr.  W.R.  Briley 

Scientific  Research  Associates,  Inc. 
P.0.  Box  498 
Glastonbury,  CT  06033 

Professor  P.  Gordon 
Caljrtef  Campus 
Department  of  Mathematics 
Purdue  University 
Hamnond,  IN  46323 

Library  (MS  135) 

NASA  Lnaqley  Research  Center 
Lanqley  Station 
Hampton,  VA  23665 

Professor  A.  Shapmaon 
Chariman,  Mechanical  Engineering 
Department 

William  M.  Rice  Institute 
Box  1892 

Houston,  TX  77701 

Technical  Library 
Naval  Ordnance  Station 
Indian  Head,  MD  20640 

°rofessor  D.  A.  Caughey 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Cornell  University 
Ithaca,  NY  14850 


Professor  E.L.  Resler 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Cornell  IJni  ve r i  sty 
Ithaca,  NY  14850 

Professor  S.F.Shen 
Sibley  School  of  Mechanical  and 
Aerospace  Engineering 
Cornell  University 
Ithaca,  NY  14350 

Li  brary 

Midwest  Research  Institute 
425  Volker  Boulevard 
Kansas  City,  MO  64110 

Dr.  J.J.  Riley 
Flow  Research,  Inc. 

P.  0.  Rox  5040 
Kent,  WA  98031 

Dr.  S.A.  Orszag 

Cambridge  Hydrodynamics,  Inc. 

54  Baskin  Road 
Lexington,  MA  02173 

Dr.  P.  Bradshaw 

Imperial  College  of  Science  and 
Technol ogy 

Department  of  Aeronaitics 
Prince  Consort  Road 
London  SW7  28 Y,  England 

Professor  T.  Cebeci 
Mechanical  Engineering  Department 
California  State  Universitv,  Long 
Beach 

Long  Beach,  CA  90340 
Dr.  H.  K.  Cheng 

University  of  Southern  California 
Department  of  Aerospace  Engineering 
University  Park 
Los  Angeles,  CA  90007 

Professor  J.  D.  Cole 
Mechanics  and  Structures  Department 
School  of  Engineering  and  Applied 
Science 

University  of  California 
Los  Angeles,  CA  90024 


Engineering  Library 

Professor  A.  Jameson 

University  of  Southern  California 

Courant  Institute  of  Mathematical 

'tox  7  7  9  2 '■? 

Sciences 

I.o s  Anqeies,  CA  9000  7 

New  York  University 

251  Mercer  Street 

Hr.  C.-M.  Ho 

New  York,  NY  10012 

Department  of  Aerospace  Engineering 

University  of  Southern  California 

Professor  G.  Miller 

'In  i  vers  i  t.y  Park 

Department  of  Applied  Science 

*  ”  * 

Los  Angeles,  CA  90007 

New  York  University 

26-36  Stuyvesant  Street 

ft  *  • 

Command  Officer 

New  York,  NY  10003 

Naval  Ordnance  Station 

Louisville,  KY  40214 

Office  of  Naval  Research 

New  Yon<  Area  Office 

Mr.  B. H.  Little,  Jr. 

715  Broadway  -  5th  Floor 

-  s 

Lockheed-Georqi a  Company 

New  York,  NY  10003 

Department  72-74,  Zone  369 

Marietta,  GA  30061 

Dr.  A.  Vagi i o-Lauri n 

LJ 

■  .* 

Professor  F.R.G.  Eckert 

Department  of  Applied  Science 

New  York  University 

University  of  Minnesota 

26-36  Stuyvesant  Street 

«/ 

241  Mechanical  Engi needing  Building 

New  York,  NY  10003 

*  X-* 

Minneapolis,  MN  66465 

Mr.  0.  Farmer 

Dr.  Gary  Chapman 

Naval  Ocean  Research  and  Development 

Mail  Stop  227-4 

Activity 

Ames  Research  Center 

Code  332 

>V 

Moffett  Field,  CA  94035 

NS TL  Station,  MS  39522 

Library 

Librarian,  Aeronautical  Library 

Naval  °ostgraduate  School 

National  Research  Council 

Monterey,  CA  93940 

Montreal  Road 

Dr.  J.  L.  Steqer 

Ottawa  7,  Canada 

rl ow  Simulations,  Inc. 

Lockheed  Missiles  and  Space  Company 

v‘* 

/3 6  Alice  Avenue 

Technical  Information  Center 

Mountain  View,  CA  94041 

3251  Hanover  Street 

.  - 

Palo  Alto,  CA  94304 

\ 

Or.  S.S.  Stahara 

"  - 

Nielsen  Engineering  arid  Research 

Commanding  Officer 

Inc. 

Office  of  Naval  Research 

510  Clyde  Avenue 

Western  Regional  Office 

-m  — 

Mountain  View,  CA  94043 

1030  East  Green  Street 

Pasadena,  CA  91106 

,N 

Engineering  Societies  Library 

s’ 

346  East  47th  Street 

Engineering  Division 

New  York,  New  York  10017 

California  Institute  of  Technology 

Pasadena,  CA  91109 

M 

w 

**  * 

U.S.  Army  Research  Office 
P.0.  Bo<  \?.?.U 

Research  Triangle,  NC  2/709 


Ufl 


> 


*.v 


'.V 


W 


3 


US 


❖ 


Li brary 

Jet  Propulsion  Laboratory 
4000  Oak  Grove  Drive 
Pasadena,  DA  91103 

Professor  H.  Liepmann 
Department  of  Aeronautics 
California  Institute  of  Technology 
Paso-iena,  CA  91109 

Mr.  L.  I.  Chasen,  MGR-MSO  Lib. 
General  electric  Company 
Missile  and  Space  Division 
P.0.  Box  BBSS 
Phi ladelphi a,  PA  19101 

Technical  Library 
Naval  Missile  Center 
°oiet  Mugu,  CA  93042 

Professor  S.  Bogonoff 
Gas  Dynamics  Laboratory 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton  University 
Princeton,  NJ  08540 

Professor  S. I.  Cheng 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton  University 
Princeton,  NJ  08540 

Dr.  J.  F.  Yates 

Aeronautical  Research  Associates 
of  Princeton,  Inc. 

50  Washington  Road 
Princeton,  NJ  03540 

Professor  L.  Sirovich 
Division  of  Applied  Mathematics 
Brown  University 
Providence,  Rl  02912 

Redstone  Scientific  Information 
Center 

Chief,  Document  Section 
Army  Missi le  Command 
Redstone  Arsenal,  AL  35309 


Editor,  Applied  Mechanics  Review 
Southwest  Research  Institute 
3500  Culebra  Road 
San  Antonio,  TX  73228 

Library  and  Information  Services 
General  Dynamics  -  CQNVAIR 
P.  0.  Box  1  i 2 3 
San  Dieqo,  CA  92112 

Dr.  R.  Magnus 

General  Dynamics  -  CONVAIR 
Kearny  Mesa  °lant 
P.0.  Box  30347 
San  Diego,  CA  92133 

Office  of  Naval  Research 
San  Francisco  Area  Office 
One  Hallidie  Plaza,  Suite  601 
San  Francisco,  CA  94102 

Li  brary 

The  RAND  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90401 

Dr.  P.  E.  Rubbert 
Boeing  Aerospace  Company 
Boeing  Military  Airplane  Development 
Qrgani zation 
P.0.  Box  3707 
Seattle,  WA  93124 

Dr.  H.  Yoshihara 
Boeing  Aerospace  Company 
P.0.  Box  3999 
Mail  Stop  41-18 
Seattle,  WA  93124 

Li  brarian 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Professor  J.H.  Ferziger 
Department  of  Mechanical  Engineering 
Stanford  University 
Stanford,  CA  943)5 


Professor  1C.  Karamcheti 
Department  of  Aeronautics  and 
Astronautics 
Stanford  University 
Stanford,  CA  94305 

9rofe$sor  0.  Runemann 
Institute  for  Plasma  Research 
Stanford  University 
Stanford,  CA  94305 

Engineering  Library 
McDonnell  Douqlas  Corporation 
Department  218,  Ruildinq  101 
P.0.  Rox  516 
St.  Louis,  MO  63166 

Or.  R.J.  Haxkinen 

McDonnell  Douqlas  Corporation 

Department  222 

P.  0.  Rox  516 

St.  Louis,  M0  63156 

Dr.  N.  Mai  noth 

Rockwell  International  Science 
Center 

1049  Camino  Dos  Rios 
P.0.  Rox  1085 
Thousand  Oaks,  CA  91860 

Li  bra  ry 

Institute  of  Aerospace  Studies 
University  of  Toronto 
Toronto  5,  Canada 

Professor  W.R.  Sears 

Aerospace  and  Mechanical  Engineering 

University  of  Arizona 

Tucson,  AZ  85721 

Professor  A.R.  Seebass 
Department  of  Aerospace  and 
Mechanical  Engineering 
University  of  Arizona 
Tucson,  AZ  85721 

Dr.  K.  T.  Yen 
Code  3015 

Naval  Air  Development  Center 
Warminster,  PA  18974 


Air  Force  Office  of  Scientific 
Research  (SREM) 

Ruildinq  1410,  Rolling  AFR 
Washington,  DC  20332 

Chief  of  Research  and  Development 
Office  of  Chief  of  Staff 
Department  of  the  Army 
Washington,  DC  20310 

Library  of  Congress 

Science  and  Technology  Division 

Washington,  DC  20540 

Director  of  Researco  (Code  RR) 
National  Aeronautics  and  Space 
Administration 
600  Independence  Avenue,  SW 
Washington,  DC  20546 

Li  brary 

National  Bureau  of  Standards 
Washington,  DC  20234 

National  Science  Foundation 
Engineering  Deivision 
1800  G  Street,  NW 
Washington,  DC  20550 

AIR  3200 

Naval  Air  Systems  Command 
Washington,  DC  20351 

AIR  950D 

Naval  Air  Systems  Command 
Washington,  DC  20361 

Code  2627 

Naval  Research  Laboratory 
Washington,  DC  20375 

SEA  03512 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

SEA  09G3 

Naval  Sea  Systems  Command 
Washington,  DC  20362 


"■  ■? ■ 


.V-i  . 


R 


>S  >fc 


■• .  -i 


V. . 


'‘-.v 

■V 

^■j 


*<» 

.■_v 


c!s 


& 


>S 

>\sV  .\V 


A  V 


m  h 

1 


8 


"H 


Or.  Charles  Watkins 
Head,  Mechanical  Engineering  Dept. 
Howard  University 
Washington,  DC  20059 

Dr.  A.L.  Slafkosky 
Scientific  Advisor 
Commandant  of  the  Marine  Corps 
(COde  AX) 

Washington,  D.  C.  20380 
Di  rector 

Weapons  Systems  Evaluation  Group 
Washington,  DC  20350 

Research  Library 
AVC0  Corporation 
Missile  Systems  Division 
201  Lowell  Street 
Wilmington,  MA  01887 


Administrative  Contracting  Officer 
Bethpage,  L. I. ,  New  York  11714 


AFAPL  ( APRC ) 

AB 

Wright  Patterson  Ai 
OHIO  45433 


Force  Base, 


Dr.  Donald  J.  Harney 
AFFDL/FX 

Wright  Patterson  Air  Force  Base, 

OHIO  45433 

FI  ow  Research 
1320  Fenwick  Lane 
Suite  401 

Silver  Spring,  MD  20910 

Scientific  Officer 
Program  Director  Fluid  Dynamics 
Mathematical  and  Information 
Sciences 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 

Attn:  Dr.  Robert  E.  Whitehead 
Code  438 

Ref:  N00014-81-C-0549 

Director,  Naval  Research  Laboratory 
Attn:  Code  2627 

Washington,  DC  20375  6  copies 


A\.\vs^ 


mL 


-vs 


"*  “W  * 


-  _  »  j.  ^  . 


